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ABSTRACT 


Three  lower  confidence  interval  estimation  procedures  for  system  reliability  of  co¬ 
herent  systems  with  cyclic  components  arc  developed  and  their  accuracy  evaluated  by 
Monte  Carlo  methods.  Each  method  uses  estimates  of  the  ratios  of  component  unreli¬ 
abilities  and  the  Poisson  approximation  to  the  binomial  distribution  to  obtain  the 
equation  for  the  lower  confidence  limit.  This  is  an  extension  of  a  method  previously 
reported  in  the  literature  which  has  been  shown  to  be  fairly  robust.  The  procedures  de¬ 
veloped  here  can  be  combined  with  similar  procedures  already  developed  for  systems 
with  continuous  components.  The  combined  procedure  may  yield  a  reasonably  accurate 
lower  confidence  interval  procedure  for  the  reliability  of  coherent  systems  with  mixtures 
of  continuous  and  cyclic  components. 


iii 


Acoeasion  F 

RTIS  GRAJbl 
DTIC  TAB 
Unannounced 
Justification 


Rv 

Distribution/ 

Availability  Codes 

•V 

lAvall  and/or 

Dlst 

i 

Speo 

lal 

(H 

□  □ 


THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  devc'.oncd  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made,  within 
the  time  available,  to  ensure  that  the  programs  arc  free  of  computational  and  logic  er¬ 
rors,  they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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Figure  1.  Block  Diagram  of  a  Wheatstone  Bridge  Structure. 


I.  BACKGROUND 


A  coherent  system  is  any  system  whose  reliability  is  not  reduced  when  the  reliability 
of  any  of  its  components  is  increased.  This  thesis  deals  with  coherent  systems  that  have 
cyclic  components.  Cyclic  components  perform  one  or  more  repetitions  of  a  single 
function  during  their  mission.  Examples  are  components  that  have  power-tum-ons  and 
power-turn-ofls,  switches  and  hydraulic  units.  Any  component  that  performs  a  function 
which  is  regarded  as  a  success  or  failure  is  a  cyclic  component.  Consequently  a  contin¬ 
uously  operating  electronic  component  becomes  a  cyclic  component  if  the  only  data  for 
that  component  is  the  number  of  successes  and  failures  in  repeated  trials.  The  number 
of  trials  to  first  failure  of  cyclic  components  have  discrete  probability  distributions. 
Numerous  approximate  interval  estimation  procedures  have  been  developed  for  system 
reliability  of  coherent  systems  with  cyclic  components.  Such  a  method  developed  by 
Mvhrc  and  Saunders  [Ref.  1:  p.37]  uses  a  likelihood  ratio  method.  C.  R.  Rao  [Ref.  2  ] 
developed  a  maximum  likelihood  method  and  Easterling  [Ref.  3  ]  proposed  a  modified 
maximum  likelihood  method.  For  special  structures,  such  as  series  systems,  exact 
methods  were  developed  by  Winterbottom  [Ref.  4:  pp.782-7S7].  Other  approximate  and 
asymptotic  methods  were  developed  by  Mann  and  Grubbs  [Ref.  5  ].  The  accuracy  of 
some  of  these  approximate  procedures  has  been  studied  for  specific  sample  size  cases  for 
structures  of  order  two  or  three  by  Mann  and  others  [Ref.  6]  and  Winterbottom  [Ref. 
4:  pp.782-787].  Woods  and  Borsting  [Ref.  7  ]  developed  an  approximate  procedure 
which  was  modified  extensively  by  Mann  and  Grubbs  [Ref.  5:  pp.335-347].  None  of 
these  discrete  interval  estimation  procedures  cited  above  can  be  readily  used  in  con¬ 
junction  with  continuous  component  data  to  obtain  interval  estimates  for  the  reliability 
of  complex  systems  that  have  mixtures  of  cyclic  and  continuous  components. 

In  this  thesis  we  attempt  to  establish  feasible  interval  estimation  methods  for  the 
reliability  of  coherent  systems  with  cyclic  components.  These  methods  have  a  common 
feature  that  allow  them  to  be  combined  with  similar  methods  that  use  continuous  data. 
The  combination  of  these  methods  may  provide  interval  estimates  for  the  reliability  of 
systems  with  cyclic  and  continuous  components.  This  method  is  an  extension  of  a 
method  developed  by  Lomnicki  [Ref.  8  ]  and  extended  by  Myhrc  and  others  [Ref.  9: 
p.2I3]l  They  compare  the  accuracy  of  their  procedure  with  the  accuracy  of  three  of  the 
procedures  cited  above  and  provide  bounds  for  error  between  the  true  lower  confidence 
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bound  and  their  approximate  lower  bound.  These  bounds  are  used  to  show  that  their 
procedure  is  fairly  robust  against  errors  in  assumptions  about  the  ratios  between  com¬ 
ponent  unreliabilities  which  are  needed  to  use  their  procedure.  The  procedure  developed 
in  this  thesis  attempts  to  capitalize  on  this  robustness  to  extend  their  method  by  esti¬ 
mating  these  ratios  from  the  test  data  .  Specifically,  let  denote  the  reliability  of  com¬ 
ponent  i  and  denote  component  unreliability  by  q,  =  1  -  Rr  Suppose  the  component 
unreliabilities  of  a  coherent  system  of  size  k  are  denoted  by  qu  q2, ... ,  qk .  Let 

qm  =  ma\{q^q2,...,qk}  (1.1) 

Qi 

and  a,  =  Then  the  reliability  function  h(<7„  q2, ... ,  qk),  of  the  system  can  be  written 
as  h (qm,  au  a2l ... ,  ak).  Since  the  system  is  coherent,  a  lower  confidence  limit  on  system 

a 

reliability  Rsx  is  given  by 

...  ~  U-2) 

* 

where  qu  is  the  corresponding  upper  confidence  limit  on  qm.  In  this  thesis  the  a,  arc  es¬ 
timated  from  the  data.  The  method  developed  by  Myhre  and  others  [Ref.  9:  pp.216-223] 
assume  the  a ,  are  known. 

Myhre  and  others  [Ref.  9:  p.223]  also  apply  this  failure  rate  ratio  concept  to  systems 
whose  components  have  exponentially  distributed  failure  times  with  failure  rates  In 

this  case, 

;.m  =  max{;.„;.2,  (i.3) 

««— T-  U-4) 

Am 

A 

and  h(yiu,  ax, ... ,  aA)  is  the  lower  confidence  limit  for  the  reliability  function 

A 

h(9.m,  au  ... ,  ak)  where  is  the  upper  confidence  limit  for  .  In  his  thesis,  Lee  [Ref. 

A 

10:  p.7]  has  demonstrated  that  if  the  same  life  data  used  to  obtain  /„  is  used  to  construct 
estimates  a,  for  the  a„  in  equation  (1.4),  then  the  approximate  lower  confidence  limit 

h(L«i,«2 . °k)  (1-5) 


is  an  accurate  approximate  interval  estimation  procedure. 


A  primary  reason  for  investigating  the  accuracy  of  this  failure  rate  ratio  method  for 
cyclic  systems  is  that  it  can  be  easily  extended  to  systems  that  have  mixtures  of  cyclic 
and  continuous  components.  If  the  procedure  investigated  here  for  systems  of  cyclic 
components  is  found  to  be  accurate,  then  it  can  readily  be  combined  with  the  procedures 
established  by  Lee  (Ref.  10:  pp.3-23]  to  obtain  an  accurate  interval  estimation  procedure 
for  the  reliability  of  a  coherent  system  with  a  mixture  of  cyclic  and  continuous  compo¬ 
nents. 

An  additional  primary  purpose  of  this  thesis  is  to  develop  an  extensive  computer 
simulation  program  that  will  provide  a  means  for  evaluating  the  accuracy  of  proposed 
interval  estimation  method*  for  the  reliability  of  coherent  systems  with  cyclic  compo¬ 
nents. 


-00 
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II.  INTERVAL  ESTIMATION  PROCEDURES  FOR  SYSTEM 

RELIABILITY 

A.  GENERAL  CONCEPTS 

Consider  a  coherent  system  with  k  cyclic  components.  Mission  tests  for  component 
/  consist  of  independent  Bernoulli  trials  administered  under  a  mission  operating  envi¬ 
ronment  with  probability  q,  of  failure  on  each  test.  Then  F„  the  number  of  failures  in  n, 
independent  mission  tests,  has  a  binomial  distribution  denoted  by  BIN {ntq).  We  shall 
refer  to  q,  as  the  unreliability  of  component  i . 

Utilizing  the  concept  first  suggested  by  Lomnicki  [Ref.  8:  p.  109  ]  and  expanded  by 
My  lire  and  others  [Ref.  9:  p.213],  the  unreliability  of  each  component  can  be  written  as 
a  fraction  of  the  unreliability  of  the  least  reliable  component;  that  is 

<h  =  aiqn  (2.1) 

where 

<lm  —  niax{^1.  i=  1,2, ,  k  and  0  <,  a,  <  1.  WeTshall  assume  that  the  q,  are  suf¬ 
ficiently  small  so  that  the  probability  distribution  of  F,  can  be  approximated  by  the 
Poisson  distribution  with  mean  n,  q,.  We  denote  this  by  the  expression 

Fi^P(n,q^.  (2.2) 

Therefore,  the  distribution  of  the  total  number  of  failures,  F  =  XT’,,  can  be  approximated 
by  the  Poisson  distribution  with  mean  q,\  that  is, 


Applying  equation  (2.1)  to  equation  (2.3)  we  obtain 


Let  h(<7„<7„ ...  ,qk)  denote  the  equation  for  system  reliability,  Rs.  Then  from 
equation  (2.1). 
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^5  ^2»  *  ®Jk)* 


(2.5) 


Since  the  system  is  coherent,  h  is  non-increasing  in  qm.  Therefore,  a  lower 
100(l-a)%  confidence  limit,  RS,LW  for  Rs  is 

'  R-S,L(3)  =  M9m.u(a)'  al*  —  *  ak)  (2-6) 


where  q„Mt)  is  any  100(1  -a)%  upper  confidence  limit  for  qm.  To  obtain  qmMt)  we  use  the 
following  well  known  result  for  the  upper  confidence  limit  of  the  mean  of  a  Poisson 
distribution,  see  [Ref.  11:  p.218].  If  X~  P{).)  then  the  upper  I00(l-a)%  confidence  limit 
for  X  using  one  observation  on  X  is 


y.a.iO+X) 
*■<*>  ~  2 


(2-7) 


Applying  this  result  to  equation  (2.4)  we  obtain  the  upper  100(1  -a)%  confidence  limit 

*  i 

for  a,  ^vhich  is  given  by 

i=i 


(2-8) 


The  corresponding  upper  confidence  limit  for  qm  is 

**.2(1  +F) 
QmjuW  ~  k 

2  a‘ 
(=1 


(2.9) 


Equation  (2.9)  is  the  expression  used  by  Myhrc  and  others  [Ref.  9  :  p.214  ]. 

The  components  of  the  vector  a  —  (a„  a2, ... ,  a*)  as  defined  by  Myhre  and  others 

[Ref.  9:  p.2 14]  are  computed  in  relation  to  the  largest  unreliability  qm.  However  any 

component  could  be  used  as  the  base  component  to  form  a  vector  of  ratios.  Choosing 

the  one  with  largest  unreliability  is  a  convenient  way  to  select  the  index  m.  In  this  thesis 

wc  will  assume  the  q,  arc  unknown.  We  will  estimate  them  from  the  data.  We  also  use 

the  data  to  determine  the  base  component,  which  we  label  component  m,  that  will  be 

used  to  form  the  ratio  estimates  a,  .  Specifically  if  A',  F,  are  the  number  of  trials  and 

/* 

resultant  failures  for  component  i,  then  q,  =  ~r  is  the  estimate  of  unreliability  q,.  I  lcrc- 
after  in  this  thesis,  we  define  the  index  m  by  the  equation 
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<Jm  ~  max(^)f  q2 . qk). 


(2.10) 


Wc  shall  define  a ,  by 


A 


/=  1,2, ...  ,k. 


(2.11) 


Note  that  qm  is  not  intended  to  be  an  estimate  of  max(^„ ... ,  qk).  It  is  nothing  more  than 
the  largest  of  the  observed  estimates  of  the  component  unreliabilities  and  provides  us 
with  the  base-component  which  we  denote  by  m  .  Hereafter,  in  this  thesis,  qm  will  denote 
the  unreliability  of  the  component  m  which  has  been  determined  by  the  definition  of  qm. 
With  the  index  m  determined  and  the  quantities  alt  a2 , ... ,  ak  computed,  we  can  com¬ 
pute  an  approximate  upper  confidence  limit  qm,u(t)  for  qm  that  corresponds  to  qmMt)  in 
equation  (2.9).  Specifically, 


A  *g,2(1  +F) 

QmMfr)  ~  k 

2  Y",  at 

/=  i 


(2.12) 


The  corresponding  approximate  lower  100(l-a)%  confidence  limit  RSLM  is 

h  (W*, . «*>•  (2-13) 


We  have  special  problems  with  equation  (2.12)  when  none  of  the  components  fail.  In 
this  case  q,  =  0,  (/'  =  1 , 2, ...  ,k)qm  =  0  and  all  a ,  are  undefined.  If  at  least  one  compo¬ 
nent  has -at  least  one  failure,  equation  (2.12)  is  well  defined.  Consequently,  equations 
(2.12)  and  (2.13)  must  be  modified  or  supplemented  to  account  for  cases  when  none  or 
perhaps  only  one  of  the  components  fail.  In  fact  we  may  want  to  modify  the  definition 
of  a,  whenever  a,  —  0;  i.e.  when  F,  =  0.  Anytime  a,  =  0,  the  N,  tests  for  component  / 
make  no  contribution  toward  the  evaluation  of  qm,uM  in  equation  (2.12).  If  information 
about  the  range  of  a,  =  •—  is  known,  it  could  be  used  to  redefine  a,  so  that  a ,  n,  in  the 
denominator  of  equation  (2.12)  would  not  be  zero.  This  is  mathematically  equivalent 
to  adding  a,  n,  tests  to  the  Nm  tests  for  component  m  without  adding  any  failures  to  Fm. 
This  is  the  one  advantage-of  the  ratio  method  and  it*is  reflected  in  both  equations  (2.12)1(W, 
and  (2.13).  It  is  this  specific  property  of  the  ratio  method  that  makes  it  particularly 
appealing  when  the  total  number  of  failures  is  small. 
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The  primary  purpose  of  this  thesis  is  to  initiate  the  development  of  alternative 
methods  for  constructing  approximate  lower  confidence  limits  RSiU))  for  the  reliability 
Rs  of  coherent  systems  that  account  for  zero  or  nearly  zero  failures  and  to  construct  a 
computer  program  that  can  be  used  to  evaluate  the  accuracy  of  these  confidence  limit 
procedures.  It  is  highly  desirable  that  these  methods  make  strong  use  of  the  ratio 
method,  whe&  several  typesjof  components  experience  failure,  because  doing  so  should 
allow  us  to  extend  these  confidence  limit  procedures  to  systems  that  have  both  cyclic  and 
continuously  operating  types  of  components. 

It  is  highly  unlikely  that  any  one  confidence  limit  procedure  will  be  reasonably  ac¬ 
curate  for  all  system  configurations;  i.e.,  for  all  reliability  functions,  h (q„ ... ,  qk).  Con¬ 
sequently,  we  have  developed  three  similar  procedures  which  are  labeled  Procedure  1,  2, 
and  3.  Each  of  these  procedures  will  be  evaluated  for  accuracy  when  applied  to  two 
radically  different  types  of  system  configurations;  namely  series  system  and  Wheatstone 
bridge  systems.  Series  systems  have  no  redundant  components.  Wheatstone  bridges  are 
highly  redundant  systems.  These  accuracy  evaluations  should  help  establish  some  pre¬ 
liminary  boundary  constraints  on  the  application  of  these  procedures  and  provide  in¬ 
sights  for  modifications  that  may  yield  more  accurate  procedures. 

The  first  evaluations  will  be  performed  for  a  series  system  of  k  independent  compo- 

*  k 

nents.  In  this  case  the  system  reliability  Rs  is  given  by  Rs~  flO  —<7,)*®  )TI(1  —  a,qm) 

i~\  /= 1 

where  q,  denotes  the  unreliability  of  component  / ,  and  m  will  be  the  index  established 
by  the  data  as  previously  discussed,  and  a,  =  .  One  expression  that  can  be  used  to 

construct  a  lower  100(l-a)%  confidence  limit  Rs,nt)  for  Rs  using  qm<uM  given  in  equation 
(2.12)  is 


k 

(=i 

In  the  following  descriptions,  and  F,  denote  the  number  of  tests  (Bernoulli  trials) 
and  the  number  of  failures,  respectively,  for  component  / ,  i  —  1 ,  2, ... ,  k. 

b.  PROCEDURE  1 

A 

This  procedure  has  three  expressions  for  Rs,u»  which  depend  on  the  test  results 
F„  F2, ... ,  Fk.  If  no  components  fail,  all  q,  =  0,  and  all  a,  arc  undefined.  This  precludes 

A 

the  use  of  equation  (2.12)  to  construct  Rs,l^  f°r  this  case.  When  no  components  fail, 
we  define 
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n*  =  min(/tt,  n2, ... ,  nk) 


(2.15) 


and  interpret  the  data  as  rr  successful  system  tests.  Rather  than  use  the  standard  ap¬ 
proximate  binomial  confidence  limit  (a)^,  we  choose  to  define  Rs,w  by 


Rs.  L(a) 


(2.16) 


This  expression  reflects  the  Poisson  approximation  to  the  Binomial  distribution  and  will 
be  slightly  conservative  when  all  F,  =  0,  but  the  other  components  to  Procedure  1  are 
suspected  to  be  slightly  optimistic.  So  the  entire  procedure  may  be  nearly  exact.  If  the 
entire  procedure  is  conservative,  we  can  always  change  equation  (2.16). 

The  second  expression  in  Procedure  1  addresses  the  case  when  all  but  one  compo¬ 
nent,  say  component  1,  have  zero  failures.  For  a  series  system,  this  would  correspond 
to  at  most  one  system  failure  out  of  nm  system  tests.  This  would  not  be  the  case  in  other 
systems  if  the  failed  component  is  redundant  in  the  system.  Here  again  we  choose 
slightly  conservative  procedure  for  the  lower  confidence  limit;  namely,  we  define  qUuM 
by 


a  Xa.2(l+F,) 

IlMi*)  2/1, 


(2-17) 


A 

We  use  equation  (2.13)  to  define  Rs,lm  for  a  general  coherent  system.  If  the  system  is 
a  series  system  we  use  equation  (2.14). 

The  third  expression  in  Procedure  1  addresses  the  case  when  at  least  two  different 
types  of  components  have  failures;  that  is,  F,  ¥=  0  for  at  least  two  values  of  /  , 


/=  1,  2, ... ,  k..  In  this  case  we  use  equation  (2.12)  to  define  qm,y{t)  where  a,  —  , 

^  h  A  A 

<7,  =  and  q„  =  max(<7„  ... ,  qk).  The  lower  confidence  limit  for  the  system  is  given  by 
equation  (2.13)  for  general  coherent  systems  or  equation  (2.14)  for  a  scries  system. 

C.  PROCEDURE  2 

Procedure  2  differs  from  Procedure  1  in  two  respects.  First,  we  assume  that  at  least 
two  different  components  of  the  system  have  at  least  one  failure.  Operationally  this 
means  that  Procedure  2  would  only  be  applied  to  data  for  which  at  least  two  types  of 
system  components  experienced  one  or  more  failures.  In  this  procedure  we  always  use 


-m 


-rtWV 
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equation  (2.1)  to  obtain  qmMm)  and  equation  (2.13)  or  (2.14)  to  compute  RSt'm  for  general 

A 

coherent  systems  or  series  systems  respectively.  Again  a ,  =  but  we  will  define  q, 

differently  for  components  with  zero  failures.  Specifically,  we  apply  a  common  scaling 
factor  K  to  increase  the  sample  size  from  n ,  to  Kn,  for  all  system  components  i  which 
experienced  zero  failures  and  record  1  failure  for  each  of  those  components.  For  ex¬ 
ample,  if  «„  and  tests  were  performed  on  component  types  1,  2  and  3  and  no  failures 
occurred  on  any  of  these  tests,  the  data  is  modified  to  show  F,  =  1 ,  F2  =  1  ,  F3  =  1  and 
new  sample  sizes  Knu  Kn2  and  Kny  The  scaling  factor  is  an  input  parameter  to  the 
computer  programs  used  to  evaluate  this  procedure;  ie.,  it  is  assumed  to  be  known.  It 
is  determined  by 


_  max(?1  ’  •••  »  ?/<)  _  ffmax 

min(9],  ...  ,  qk)  ^min 


(2.18) 


where  the  q,  are  the  component  unreliabilities.  Operationally,  this  would  mean  that  the 
user  of  this  procedure  must  determine  an  estimate  for  this  ratio.  Previous  data  is  often 
available  to  provide  this  estimate.  Such  estimates  can  sometimes  be  constructed  from 
Department  of  Defense  (DOD)  documents  that  provide  a  variety  of  quality  data  char¬ 
acteristics  for  hardware  purchased  in  accordance  with  prescribed  standards  and  specifi¬ 
cations;  c.g.,  MILHDBK-2I7E  and  its  referenced  documents.  Also  the  reliability  values 
for  components  used  in  reliability  prediction  analysis  are  usually  derived  from  some  of¬ 
ficial  source  which  could  bff'satisfactbry  for  such  estThiates.  A  reliability  prediction 
analysis  is  usually  required  in  any  major  DOD  system  acquisition  program.  If  a  satis¬ 
factory  estimate  for  K  cannot  be  obtained,  then  Procedure  2  should  not  be  used.  Of 
course  one  of  the  purposes  of  this  thesis  is  to  investigate  the  accuracy  of  this  procedure 
when  the  scaling  factor  K  is  chosen  “correctly";  i.e.,  according  to  equation  (2.18). 
Thinking  of  qmtx  and  qm[n  as  failure  rates,  —  and  —  denote  the  expected  number  of 
tests  to  first  failure  for  their  respective  components.  Then  the  value  K  in  equation  (2.18) 
is  the  ratio  of  these  two  expected  number  of  tests  to  obtain  one  failure  for  each  compo¬ 
nent.  That  is, 
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F or  each  system  component  /  that  experiences  zero  failures,  the  corresponding  esti¬ 
mate  q,  for  <7,  is  given  by 


1 

Kn(  ‘ 


Then 


3l 

A 

<7m 


and 


2 

a  Xa,2(l+f) 

i=0 


(2.20) 


(2.21) 


(2.22) 


where  F  is  the  number  of  actual  failures  observed  before  the  data  is  adjusted.  Thus  the 

data  is  adjusted  to  obtain  estimated  values  a,  different  from  zero  for  those  components 

* 

that  have  no  failures.  The  value  for  is  given  by  equation  (2.14).  Note  that  if  all 

A 

k  types  of  components  have  at  least  one  failure,  the  resulting  value  of  RStU is  the  same 

■■V'-  "■***'  ■*>*> 

under  Procedures  1  and  2. 

D.  PROCEDURE  3 

Procedure  3  differs  from  Procedure  1  in  only  one  respect.  We  scale  the  sample  sizes 
as  we  did  in  Procedure  2.  That  is,  for  each  system  component  i  that  has  no  failures,  we 
estimate  q,  from  equation  (2.20).  Equations  for  a„  qmM  and  Rs.l>,)  arc  given  by  equations 
(2.21),  (2.22)  and  (2.13)  or  (2.14)  respectively. 

Procedure  3  does  not  require  that  at  least  one  failure  have  occurred  on  at  least  two 
different  types  of  components  as  required  in  Procedure  2.  Operationally,  Procedure  3 
can  be  used  for  all  sets  of  data  including  those  sets  where  no  failures  occur. 

E.  BRIDGE  SYSTEM 

We  define  a  Wheatstone  bridge  system  by  the  reliability  block  diagram  shown  in 
Figure  1  on  page  11.  If  p,  denotes  the  reliability  of  component  i,  i=  1,2, ...  ,5,  and 
1  —  />,  =  qt  =  a ,  qm  where  qm  -  max{</|,  q2 , ... ,  q%),  then  system  reliability  Rs  is  given  by 
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1 


Figure  1.  Block  Diagram  of  a  Wheatstone  Bridge  Structure. 


Rs  ~  h(<7m>  flj, ... ,  as) 

=  1  “  <?m(<h«2  +  aAas)  -  ql{axa3a5  +  a2a3aA) 

4  \ 

+  4"  ^1^2^3^5  "b  ^]^2^4^S  "h  4"  Q2(l3O^Q^) 

“  2(]m(“\Cl2a3Ws) 

(see  Myhre  and  others  (Ref.  9:  p.215])/'  Then  the  equation  for  the  lower  100(l-ct)% 
confidence  limit  on  this  system  is 

Rs,  L(y)  =  h(9m,  «(,)•  a\ »  ...  *  fls)  (2-24) 

where  q„Mt)  and  a,  arc  defined  as  in  Procedures  1,  2  or  3.  Only  Procedure  1  was  evalu¬ 
ated  for  the  bridge  system  and  reported  in  this  thesis. 
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III.  SIMULATION  PROCEDURE 


Standard  computer  simulation  methods  were  used  to  evaluate  the  accuracy  of  the 
three  interval  estimation  methods  discussed  in  section  II.  The  input  parameters  needed 
to  run  the  computer  programs  for  each  of  the  three  interval  estimation  procedures  are 
as  follows: 

•  k:  number  of  components  in  the  system 

•  H~  («n  ••• »  «*)•  vector  of  sample  sizcs(mission  tests) 

•  H  -  (?i»  ft*  » <?*):  vector  of  component  unreliabilities 

•  a:  level  of  confidence 

The  NON-IMSL  uniform  random  number  generator  SRND  was  used  to  simulate 

A 

the  outcomes  of  mission  tests  on  the  k  components.  A  total  of  £/j,  uniform  random 
numbers  were  generated  and  transformed  into  ones  or  zeroes  as  follows:  the  ni  numbers 
„rn,  xl2, ... ,  x,„(  in  the  ;th  block  of  uniform  random  numbers  were  transformed  by  the  ex¬ 
pression 


YL,  =  1  if  Xfj  <  qt 
Yy  =  0  if  Xy  >  qt 

for  j  —  1,2, ... ,  n„  and  i  =  1,2, ... ,  k,  where  q,  is  the  input  parameter  denoting  the  unre¬ 
liability  of  component  /.  Then  F,  =  ’Z.Y,j  denotes  the  number  of  failures  in  n,  independent 
Bernoulli  trials.  This  set  of  data  is  used  to  compute  the  values  of  q„  a„  qm  ltW  and  Rs,Li%) 
for  each  of  the  three  interval  estimation  methods.  Thus  one  realization  of  the  lower 

A  A  A 

bound  estimate  ,  for  the  random  variable  Rs  L(^  is  obtained  for  each  set  of  ^n,  uni- 

/«0 

form  random  numbers  generated.  Each  realization  described  was  replicated  1000  times 

A 

in  order  to  generate  a  simulated  empirical  population  of  this  random  variable  RS  L . 
The  vector  of  replications  thus  obtained  was  sorted  from  smallest  to  largest  to  obtain 

A  A 

Rlmad'  »  ^lcwioco)-  Then  the  100(1 -a)%  approximate  lower  confidence  bound  on  the 
true  system  reliability  Rs  is  the  1000(l-a)th  element  of  the  sorted  vector  of  replications; 

A 

namely,  iooo<i-»»  • 

The  simulated  true  confidence  level  is  calculated  by  finding  the^clcmcnt  of  the  vector 
of  replications  which  is  closest  to  Rs  and  observing  its  ipdex  number  j.  Then  the  sirnu- 
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lated  true"  confidence  level  is  x  100.  If  there  is  a  sequence  of  indices  with  the 
same  value  closest  to  Rs  ,  we  chose  the  smallest  index  j.  This  may  yield  an  arbitrarily 
small  value  for  our  recorded  "true"  confidence  level. 

In  addition  to  the  external  subroutines  programmed  by  the  author  and  to  the  ones 
mentioned  above,  the  following  subroutines  were  used: 

•  1MSL  subroutine  MDCH1  was  used  to  calculate  the  value  of  a  chi-square  random 
variable  given  the  quantile  and  the  degrees  of  freedom. 

•  NON-IMSJL  subroutine  SHSORT  >yas  used  to  sor^the  arrays  of  system  reliability 
,  estimates  in  ascending  order  in  order  to  obtain  the  appropriate  order  statistic. 

•  IMSL  subroutine  USMNMX  was  used  to  extract  the  minimum  and  maximum 
values  of  the  vector  of  mission  tests  for  each  component  and  of  the  vector  of  input 
component  unreliabilities. 

Each  set  of  input  parameters  defines  a  simulation  run,  or  a  case.  The  following  ta¬ 
bles  specify  the  sets  of  parameters  used  for  each  case  that  was  simulated. 


Table  1.  INPUT  NUMBER  OF  MISSION  TRIALS  N,  AND  UNRELIABILITIES 
Q,  FOR  A  SERIES  SYSTEM.  CASES  1-3.  _  _ 


Component  no. 

«,  by  Case  Number 

Input  qi 
(cases  1-3) 

1 

2 

3 

I 

100 

50 

30 

0.0200 

2 

30 

30 

25 

0.0100 

.  3 

15 

10 

20 

0.0050 

4 

10 

10 

10 

0.0050 

5 

s  !  5 

I 

5 

0.0050 

Table  2.  INPUT  NUMBER  OF  MISSION  TRIALS  N,  AND  UNRELIABILITIES 
_  Q,  FOR  A  SERIES  SYSTEM.  CASES  4-6. 


Component  no. 

> 

n,  by  Case  Number 

Input  q , 
(cases  4-6) 

4 

5 

6 

1 

100 

50 

30 

0.0100 

2 

30 

30 

25 

0.0050 

3 

15 

10 

20 

0.0025 

4 

10 

10 

10 

0.0025 

5 

5 

5 

5 

0.0025 

Table  3.  INPUT  NUMBER  OF  MISSION  TRIALS  N,  AND  UNRELIABILITIES 
Q,  FOR  A  SERIES  SYSTEM.  CASES  7-9. 


Component  no. 

/?,  by  Case  Number 

Input  <7, 
(cases  7-9) 

7 

8 

9 

1 

100 

50 

30 

0.0200 

2 

30 

30 

25 

0.0100 

3 

15 

15 

20 

0.0050 

4 

10 

5 

15 

0.0050 

5 

10 

5 

10 

0.0050 

6 

10 

f 

kJ 

5 

0.0050 

7 

10 

5 

5 

0.0050 

8 

10 

5 

5 

0.0050 

9 

10 

5 

5 

0.0050 

10  „ 

10 

5 

0.0050 
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Table  6.  INPUT  NUMBER  OF  MISSION  TRIALS  N,  AND  UNRELIABILITIES 
_ Q,  FOR  A  WHEATSTONE  BRIDGE.  CASES  16-18. 


Component  no. 

n,  by  Case  Number 

Input  q, 
(cases  16-18) 

16 

17 

18 

1 

100 

50 

30 

0.0500 

2 

30 

30 

25 

0.0350 

3 

15 

10 

20 

0.0350 

4 

10 

10 

10 

0.0350 

5 

5 

5 

5 

0.0350 

The  values  of  the  scaling  factor  K  used  in  Procedures  2  and  3  are  as  follows: 


Table  7.  VALUES  OF  SCALING  FACTOR  A'. 


Case  Number 

K 

1-12 

4 

13-18 

1.43 

These  values  are  computed  from  the  assigned  input  parameters  for  these  cases;  c.g., 


4  = 


.01 


.0025 


from  table  4  and  1.43=*  r  from  table  6. 

.030' 
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IV.  SIMULATION  RESULTS 


The  results  of  the  simulation  runs  for  the  twelve  series  system  cases  are  presented  in 
Tables  8  through  19.  The  results  of  the  simulation  runs  for  the  six  bridge  system  cases 
are  presented  in  Tables  20  and  21.  The  accuracy  of  a  particular  interval  estimation 
procedure  for  each  case  can  be  assessed  from  these  tables.  If  the  interval  estimation 

A 

procedure  is  exact,  the  two  values  Rs  and  /?*.<.>,  icoo<i-«)  would  be  equal  and  the  True  Con¬ 
fidence  Level  value  would  be  lOO(l-a).  For  example,  consider  Table  8,  case  2  and 
a  =  .20.  True  system  reliability,  Rs  is  0.95572  and  RL(s)i  10oo(i-«>  =  0.96398.  Thus  the  sim- 

A 

ulations  indicate  that  the  80th  percentile  point  of  the  distribution  of  RL(t)  for  Procedure 
2  for  case  2  is  0.964  instead  of  0.956.  This  error  reflects  the  inaccuracy  of  the  approxi¬ 
mate  80%  lower  confidence  limit  Rs,u 2o>  as  defined  under  Procedure  2.  The  table  also 
indicates  that  Rs  s  0.95572  is  approximately  the  77.4  percentile  point  of  the  distribution 
of  Rs,hm)  instead  of  the  80th  percentile  point. 

The  number  of  components  in  the  system  that  had  zero  failures  in  each  of  the  1000 
replications  was  recorded.  The  average  of  the  1000  values  is  displayed  in  the  last  column 
of  the  table. 

■  mt'  ■(****'  -ftf* 

A.  SERIES  SYSTEM 


Tabic  8.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
2,  CASES  1-3  _  . _ 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

*5 

A 

Rut),  joooci—*) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

1 

a  =  .20 

0.95572 

0.93563 

93.20 

3.71 

a  =  .05 

0.95572 

0.939S1 

100.00 

3.71 

a  =  .20 

0.95572 

0.96398 

77.40 

3.98 

a  =  .05 

0.95572 

0.93366 

100.00 

3.98 

3 

a  =  .20 

0.95572 

0.96633 

65.40 

4.16 

a  =  .05 

0.95572 

0.93793 

100.00 

4.16 
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increase. 


Table  9.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
2,  CASES  4-6 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Ru  i).  iooo(i-«) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

4 

MBRM 1 

0.97768 

0.96737 

100.00 

4.15 

a  =  .05 

0.97768 

0.93981 

100.00 

4.15 

5 

o 

Cj 

II 

8 

0.97768 

0.96398 

100.00 

4.40 

a  =  .05 

0.97768 

0.93366 

100.00 

4.40 

6 

a  =  .20 

0.97768 

0.96633 

100.00 

4.55 

O 

II 

8 

0.97768 

0.93793 

100.00 

4.55 

Table  10.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
2,  CASES  7-9 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Ru*).  I000<1— *> 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

7 

a  =  .20 

0.93206 

0.93884 

77.30 

8.49 

a  =  .05 

0.93206 

0.94091 

94.30 

8.49 

8 

a  =  .20 

0.93206 

0.88770 

81.00 

8.86 

a  =  .05 

0.93206 

0.89162 

.100.00 

8.86 

9 

a  =  .20 

0.93206 

0.94925 

70.60 

9.01 

a  =  .05 

0.93206 

0.90729 

100.00 

9.01 

IS 


increase. 


Table  11.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
2,  CASES  10-12 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

*5 

A 

loocKi— «> 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

10 

a  =  .20 

0.96552 

0.96788 

75.30 

9.03 

a  =  .05 

0.96552 

0.94091 

100.00 

9.03 

11 

a  =  .20 

0.96552 

0.94046 

100.00 

9.36 

a  =  .05 

0.96552 

0.89162 

100.00 

9.36 

12 

R 

II 

O 

0.96552 

0.94925 

100.00 

9.48 

a  =  .05 

0.96552 

0.90729 

100.00 

9.48 

Table  12.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
1,  CASES  1-3 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Rl{i).  1000(l-«) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

1 

a  =  .20 

0.95572 

0.95721 

65.20 

3.70 

a  =  .05 

0.95572 

0.95254 

100.00 

3.70 

5 

a  =  .20 

0.95572 

0.94011 

100.00 

3.97 

a  =  .05 

0.95572 

0.9050S 

100.00 

3.97 

3 

a  =  .20 

0.95572 

0.90018 

100.00 

4.14 

a  =  .05 

0.95572 

0.841  SI 

100.00 

4.14 

19 


increase. 


Table  13.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
1,  CASES  4-6 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

*5 

A 

1000(1— «) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

4 

a  =  .20 

0.97768 

0.97005 

100.00 

4.14 

a  =  .05 

0.9776S 

0.95254 

100.00 

4.14 

5 

ot  =  .20 

0.97768 

0.94011 

100.00 

4.38 

a  =  .05 

0.97768 

0.9050S 

100.00 

4.38 

6 

a  =  .20 

0.97768 

0.88021 

100.00 

4.53 

a  =  .05 

0.97768 

0.84181 

100.00 

4.53 

Table  14.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 


1.  CASES  7-9 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Rud.  1000(l-«) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

7 

a  =  .20 

0.93206 

0.95721 

59.80 

8.48 

a  =  .05 

0.93206 

0.95254 

74.40 

8.4S 

8 

w 

II 

io 

o 

0.93206 

0.94011 

75.90 

8.85 

M 

II 

O 

Ul 

0.93206 

0.90508 

100.00 

8.85 

9 

_ 

a  =  .20 

0.93206 

0.88021 

100.00 

8.99 

a  =  .05 

0.93206 

0.84181 

100.00 

S.99 

increase. 


Table  15.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
1,  CASES  10-12 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Rut),  icooa-j) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

40 

a  —  .20 

0.96552 

0.97005 

74.40 

9.02 

«  =  .05 

0.96552 

0.95214 

100.00 

9.02 

11 

a  =  .20 

0.96552 

0.94011 

100.00 

9.35 

a  =  .05 

0.96552 

0.90508 

100.00 

9.35 

12 

a  =  .20 

0.96552 

0.88021 

100.00 

9.44 

a  =  .05 

0.96552 

0.84181 

100.00 

9.44 

Table  16.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
3,  CASES  1-3 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Rl{i),  100<XI-i) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 .... 

I 

a  =  .20 

0.95572 

0.93563 

93.20 

3.70 

a  —  .05 

0.95572 

0.939S1 

100.00 

3.70 

2 

a  =  .20 

0.95572 

0.96398 

77.30 

3.97 

a  =  .05 

0.95572 

0.93366 

100.00 

3.97 

3 

a  =  .20 

0.95572 

0.96633 

65.40 

4.14 

a  =  .05 

0.95572 

0.93793 

100.00 

4.14 

21 


increase. 


Table  17.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
3,  CASES  4-6 


No.  of 
Comprits 

Case 

no. 

Signif¬ 

icance 

Level 

,“*>  Rs 

A 

■^A(f£lOOO(I-«) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

4 

a  =  .20 

0.97768 

0.96737 

100.00 

4.14 

a  =  .05 

0.97768 

0.93981 

100.00 

4.14 

5 

a  =  .20 

0.97768 

0.96398 

100.00 

4.38 

a  =  .05 

0.97768 

0.93366 

100.00 

4.38 

6 

a  =  .20 

0.97768 

0.96633 

100.00 

4.53 

a  =  .05 

0.97768 

0.93793 

100.00 

4.53 

Table  18.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
3,  CASES  7-9 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

Rl(').  1000(1— •) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

7 

a  =  .20 

0.93206 

0.93SS4 

77.30 

8.48 

a  =  .05 

0.93206 

0.94091 

94.30 

8.4S 

8 

o 

rt 

II 

« 

0.93206 

0.88770 

80.80 

8.85 

a  =  .05 

0.93206 

0.89162 

100.00 

8.S5 

9 

a  =  .20 

0.93206 

0.94925 

70.40 

8.99 

a  =  .05 

0.93206 

0.90729 

100.00 

8.99 

increase. 


Table  19.  SUMMARY  OF  RESULTS  FOR  A  SERIES  SYSTEM-PROCEDURE 
3,  CASES  10-12 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

Rs 

A 

iooo<j-«) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

10 

10 

» 

II 

o 

0.96552 

0.96788 

75.30 

9.02 

a  =  .05 

0.96552 

0.94091 

100.00 

9.02 

11 

R 

II 

ij 

o 

0.96552 

0.94046 

100.00 

9.35 

a  =  .05 

0.96552 

0.89162 

100.00 

9.35 

12 

R 

II 

O 

0.96552 

0.94925 

100.00 

9.44 

a  =  .05 

0.96552 

0.90729 

100.00 

9.44 

B.  BRIDGE  SYSTEM 


Table  20.  SUMMARY  OF  RESULTS  FOR  A  BRIDGE  SYSTEM-PROCEDURE 
r  "  1,  CASES  13-15  _  ' 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

_ 

Rs 

A 

Run),  1000(1—®) 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

13 

R 

II 

tJ 

0.9S742 

0.99372 

56.6 

1.64 

a  =  .05 

0.98742 

1.00000 

71.50 

1.64 

1 

14 

a  =  .20 

0.9S742 

0.99332 

59.30 

1.7S 

a  =  .05 

0.9S742 

1.00000 

77.20 

1.78 

15 

a=  .20 

,.0.98742 

0-99333 

59.60 

1.57 

a  =  .05 

0.98742 

1.00000 

SO.  10 

1.57 

23 


increase. 


Table  21.  SUMMARY  OF  RESULTS  FOR  A  BRIDGE  SYSTEM-PROCEDURE 
1,  CASES  16-18 


No.  of 
Compnts 

Case 

no. 

Signif¬ 

icance 

Level 

R* 

A 

Rus),  1000(1— ») 

True  Con¬ 
fidence 
Level 

Avg.  no. 
of  comp, 
with  no 
failure 

5 

16 

a  =  .20 

0.94S16 

0.96354 

62.60 

0.77 

a  =  .05 

0.94816 

0.97426 

77.70 

0.77 

17 

a  =  .20 

0.94816 

0.96122 

64.10 

0.94 

a  =  .05 

0.94816 

0.96740 

82.80 

0.94 

18 

a  =  .20 

0.94816 

0.9585S 

69.70 

0.73 

a  =  .05 

0.94816 

0.97125 

86.30 

0.73 

V.  SUMMARY  AND  CONCLUSIONS 


The  results  in  Tables  S  through  19  in  section  IV  are  quite  mixed  with  respect  to  the 
accuracy  across  all  twelve  cases.  However,  these  results  definitely  indicate  that  any  one 
of  the  three  interval  estimation  procedures  offers  excellent  potential  for  a  reasonably 
accurate  lower  confidence  interval  estimation  method  for  system  reliability.  The  accu¬ 
racy  of  approximate  discrete  interval  estimation  methods  for  system  reliability  in  general 
depend  on  the  amount  of  testing  relative  to  the  expected  number  of  failures,  Jin,qr 
Woods  and  Borsting  [Ref.  7]  for  example,  show  their  lower  confidence  limit  procedure 
to  be  quite  accurate  for  cases  for  which  the  expected  number  of  failures  is  greater  than 
5.  The  expected  number  of  failures  for  some  of  the  cases  simulated  in  this  thesis  are  as 
follows: 


Table  22.  EXPECTED  NUMBER  OF  FAILURES 


Case  No.: 

~  1 

10  - 

►  13 

14 

16 

2>. 

2.45 

1.225 

2.725 

1.362 

14.2 

8.85 

6.8 

Consequently,  any  cases  other  than  those  numbered  13,  14,  and  16  should  not  be  ex¬ 
pected  to  be  highly  accurate  if  these  three  procedures  have  features  similar  to  those  of 
the  Woods-Borsting  procedure  [Ref.  7  :  p.  1  ].  Actually,  in  light  of  the  small  values  of 
jlrt.q,  for  cases  1,  4,  7,  and  10,  it  is  somewhat  surprising  that  the  procedures  arc  as  ac¬ 
curate  as  they  arc  at  the  80%  level  of  confidence.  The  results  for  cases  13,  14,  and  16 
arc  quite  good  at  the  80%  level  of  confidence  for  the  bridge  system.  Some  of  the  pa¬ 
rameter  sets  for  the  cases  run  were  designed  to  detect  instances  where  the  methods  might 
not  be  accurate. 

The  three  procedures  analyzed  in  this  thesis  should  be  modified  and  combined  to 
establish  one  general  procedure.  The  corresponding  range  of  parameter  sets  for  which 
the  procedure  is  accurate  should  be  established.  1  his  range  will  likely  be  a  function  of 
v„  .y  This  kind  of  a  procedure  could  be  useful  in  many  practical  settings.  Moreover, 


it  could  be  combined  with  the  method  developed  in  the  thesis  by  Lee  [Ref.  10:  pp.3-23). 
Additional  analysis,  and  simulation  need  to  be  performed  to  establish  such  a  procedure. 


APPENDIX  A.  INTERVAL  ESTIMATION  PROCEDURE  1 

PROGRAM  ZFYSCN 


* 

TITLE:  BINOMIAL  INTERVAL  ESTIMATION  PROCEDURE  * 

ZERO  FAILURES  ALLOWED;  NO  SCALING  * 

AUTHOR:  E.  F.  BELLINI,  LT,  USN  * 

DATE:  NOV  89  * 

* 


THIS  PROGRAM  COMPUTES  THE  TRUE  CONFIDENCE  LEVEL  FOR  THE  ESTIMATE  * 
RELIABILITY  OF  A  SERIES  AND  BRIDGE  SYSTEM  GIVEN  THE  RELIABILITY  * 
OF  THEIR  COMPONENTS  * 


IN  ITS  PRESENT  CONFIGURATION  THIS  PROGRAM  IS  SET  UP  TO  RUN  12  * 
TIMES  EACH  TIME  PRODUCING  1000  REPLICATIONS  USING  A  DIFFERENT  * 
SET  OF  INPUT  DATA.  RUN  THE  PROGRAM  FROM  CMS  BY  TYPING  ’B1  EXEC'.* 
THE  REXX  EXEC  PROGRAM  * 
'  B1  ’  CALLS  THE  INPUT  FILESf^TO  BE  READ  XflD  NAMES  THE  12  OUTPUT  * 
FILES  RESULTING  FROM  THE  12  CONSECUTIVE  RUNS.  BY  EDITING  THE  * 
INDEX  COUNTERS  I,  J,  K  OF  THE  ,B1’  EXEC  ONE  CAN  RUN  ANY  USER-  * 
SPECIFIC  RUN  FROM  JUST  ONE  RUN  TO  ALL  12.  * 


VARIABLES  USED 

AH AT I  : 
AI  : 
ALFA  : 
BIGF  : 
CHISQ  : 
C1C15  : 

DEGFR  : 
DELBRG  : 
DELSTR  : 
DELTAR  : 
DIFF  : 
EPS  : 
ERROR  : 
FAILS  : 
FI  : 
FLAG  : 
INC  : 
KEY1  : 
KEY  2  : 

KEY3  : 
KEY4 
KK  : 
LOOP  : 
MAXALF  : 
MAXREP  : 
MAXRUN  : 
MSTRQ  : 
MULT  : 


Vf 

ie 


WEIGHT  ESTIMATES  FOR  EACH  COMPONENT  * 

INPUT  WEIGHTS  FOR  EACH  COMPONENT  * 

LEVELS  OF  SIGNIFICANCE  * 

TOTAL  NO.  OF  FAILURES  FOR  EACH  REPLICATION  * 

CHI-SQUARE  RANDOM  VARIABLE  VALUE  * 

FORMAT  LABEL  * 

DEGREES  OF  FREEDOM 

DIFFERENCE  FOR  BRIDGE  SYSTEM  * 

DIFFERENCE  FOR  SERIES  SYSTEM-  CLOSED  FORM 
DIFFERENCE  FOR  SERIES  SYSTEM  * 

DIFFERENCE  (TRUE  REL.  -  ESTIMATED  REL. )  * 

SMALL  QUANTITYC CONSTANT)  * 

PARAMETER  FOR  IMSL  ROUTINE  *  * 

COUNTS  NO.  OF  REPLICATIONS  WITH  AT  LST.  1  FAILURE  * 
NO.  OF  FAILURES  FOR  EACH  COMPONENT (ALL  MISSION  TST)* 
1  IF  ALL  COMP.  HAVE  SAME  NO.  OF  MISSION  TESTS  * 

INCREMENT  STEP  SIZE  FOR  ROUTINE  USMNMX  * 

ARRAY  OF  INDECES  I OR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  SIZING  PARAMETER  FOR  THE  MAX  NO  OF  COMPONENTS* 
COUNTS  NO.  OF  REPLICATION  PERFORMED  * 

MAX  NO.  OF  SIGNIFICANT  LEVELS  DESIRED( ARRAY  SIZING)* 
MAX  NO.  OF  REPLICATIONS  * 

MAX  NO.  OF  FROGRAM  ITERATIONS  ALLOWED  * 

MASTER  UNRELIABILITY(USED  WITH  Al'S  TO  CALC.  Ql’S)  * 
MULTIPLIER  FOR  RANDOM  NO.  GENERATOR  SRN1  * 


k 

* 

k 

k 

k 

it 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 


k 


N  : 
NIMAX  : 
NIM1M  : 
NINDX  : 
NIREAL  : 
NMAX  : 
NPRNT  : 
PRNT  : 
QHATI  : 
QHTMAX  : 
QHTUPR  : 
QI  : 
QINDX  : 
QUANTL  : 
REPSHD  : 
"RHTSTR  : 
RS  : 
RSBRDG  : 
RSHAT  : 
RSHTBR  : 
SEED  : 
SELCTA  : 
SELCTB  : 
SORT  : 
SUMNAI  : 
TEMP  : 
TOTREP  : 
TRANBR  : 
TRANSQ  : 
TRANSR  : 
TRIALS  : 
TRNSTR  : 
TRUQNT  : 
UNIRV  : 
ZFAILS  : 
ZFPREP  : 


* 

★ 


NO.  OF  MISSION  TEST  FOR  EACH  COMPONENT 
MAX  NO.  OF  MISSION  TESTS 
MIN  NO.  OF  MISSION  TESTS 

INDEX  NO.  OF  MAX  NO.  OF  MISSION  TESTS  * 

NO.  OF  MISSION  TESTS  TRANSFORMED  TO  REAL  * 

MAX  NO.  OF  MISSION  TESTS  FOR  OUTPUT  CONTROL  * 

FLAG  FOR  DETAILED  REPORT  OUTPUT  * 

SAME  AS  ABOVE (PARAMETER)  * 

UNRELIABILITY  ESTIMATES  FOR  EACH  COMPONENT  * 

LARGEST  QHATI  * 

UPPER  LIMIT  ON  SYSTEM  UNRELIABILITY  * 

INPUT  UNRELIABILIY  FOR  EACH  COMPONENT  * 

INDEX  * 

QUANTILE  * 

REPLICATIONS  HEADING  FORMAT  NUMBER  * 

SERIES  SYSTEM  RELIABILITY  ESTIMATE( CLOSED  FORM)  * 
TRUE  SERIES  SYSTEM  RELIABILITY  * 

TRUE  BRIDGE  SYSTEM  RELIABILITY  * 

SERIES  SYSTEM  RELIABILITY  ESTIMATE  * 

BRIDGE  SYSTEM  RELIABILITY  ESTIMATE  * 

PARAMETER  * 

SIGNIFICANCE  LEVEL  SELECTION  * 

SIGNIFICANCE  LEVEL  SELECTION  * 

PARAME1ER  FOR  ROUTINE  SRND  * 

SUM  OF  THE  PRODUCT  OF  Nl'S  AND  Al'S  * 

TEMPORARY  ARRAY  * 

TOTAL  NUMBER  OF  PROGRAM  ITERATIONS  * 

TEMPORARY  ARRAY  * 

TEMPORARY  ARRAY  * 

TEMPORARY  ARRAY  * 

BERNOULLI  TRIALS  ARRAY  (2-DIM)  * 

TEMPORARY  ARRAY  * 

TRUE  QUANTILE  * 

UNIFORM  RANDOM  DEVIATES  (2-DIM)  * 

TOTAL  NUMBER  OF  REPLICATIONS  WITH  ZERO  FAILURES  * 
NO.  OF  COMPNTS.  WITH  ZERO  FAILURES  PER  REPLICATION  * 

* 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkickkkkkkkkkkkkkkkkkkkitkkkitk 


PARAMETER  (KK=10 ,MAXALF=2 ,NPRNT=0) 

PARAMETER  (MAXREP=1000 ,  MAXRUN=2000,  EPS=. OOOOOl) 

REALM  UNIRV ( 15 , 1000) ,TEMP( 1000) ,QI(KK) ,AI(KK) ,AHATI(KK) 

REALM  QHATI(KK),  NMAX,  NNMAX,  QHTMAX,  CHISQR(5,5),  ALFA(MAXALF) 
REALM  DF( 5 ) , AALFA( 5 ) , SUMNAI ,RSHAT( MAXALF , MAXREP) , RS 
REALM  KEYl(MAXRF.P)  ,KEY2(MAXREP)  ,KEY3( MAXREP)  ,TRNSTR( MAXREP) 
REALM  DEGFR( MAXREP ) ,  QHTUPR( MAXALF , MAXREP ) , CHI SQ( MAXALF , MAXREP) 
REALM  QUPA1( MAXREP) ,  QUPA2 ( MAXREP ) ,RHTSTR( MAXALF, MAXREP) 

REALM  DELTAR( MAXALF),  TRANSQ( MAXREP) , TRANSR (MAXREP) ,DIFF(MAXR£P) 
REALM  DELSTR( MAXALF) , NIMIN, NIMAX .NIREAL(KK) 

REALM  RSHTBR( MAXALF, MAXREP)  ,DELBRG( MAXALF)  ,XF.Y4( MAXREP) 

REALM  TRANBR( MAXREP ) ,  RSBRDG  ,MSTRQ 
REALM  ZFPREP 

INTEGER  SEED,  MULT,  SORtT  TRIALS ( 15  ,'f000)  ,  BIQF,  FI(.KK)  ,  N(KK) 
INTEGER  NINDX,  QINDX,  ERROR,  REPS,  SELCTA,  SELCTK ,  TOTREP 
INTEGER  C1C15 ,  REPSHD,  SELCTB,  ALF,  FLAG,  LOOP, PRNT 


2S 


INTEGER  QUANTL(MAXALF) ,  TRUQNT(MAXALF) ,ZFAILS, FAILS, INC 

DATA  SEED/123457/,  MULT/1/,  INC/1/ 

DATA  AALEA/ . 01,.  05 , .  9 , .  95 , .  99/ ,  DF/1, 5, 10,30,40/ 

DATA  ALFA/.  20,.  050/ 

DATA  S0RT/0/ 


ASSIGN  8  TO  C1C15 
ASSIGN  9  TO  REPSHD 

*  CALL  COffPRS 
*  PRNT  =  NPRNT 


DO  12  1=1, KK 

AI(I)  =  9999. 

N( I)  =  99999999 
12  CONTINUE 

READ(03,*)K,MSTRQ 

DO  11  1=1, K 
READ(03,*)  AI( I) ,N( I) 
11  CONTINUE 


IFCK.NE.5)  THEN 

WRITE( 1 , ' ( ' 'WARNING:  BRIDGE  STRUCTURE  ' ' , 

+  '  ' ONLY  USES  THE  FIRST  5  COMPONENTS")’) 

ELSE 
END  IF 

***//  INITIALIZE  THE  QHTUPR  ARRAY  OF  UNRELIABILITY  REPLICATIONS,  //*** 

*  RSHAT  ARRAY  OF  ESTIMATE  SYSTEM  RELIABILITY  REPLICATIONS  * 

*  AND  RHTSTR  ARRAY  OF  EST.  SYST.  REL.  FOR  A  SERIES  SYST  WHEN  * 

***//  ALL  THE  COMPONENT  MISSION  TESTS  ARE  EQUAL  IN  NUMBER  //*** 


173 

172 


DO  172  ALF=1 ,MAXALF 

DO  173  REPS=1 ,MAXREP 

QHTUPR ( ALF , RE  P  S )  =  0. 
RSHAT( ALF , REPS )  =  0. 
RHTSTRC ALF, REPS)  =  0. 
RSHTBRC ALF, REPS)  =  0. 
CONTINUE 
CONTINUE 


***//  SET  FLAG  TO  1  IF  ALL  COMPONENTS  HAVE  SAME  NO.  OF  MISSION  TESTS**** 


FLAG-1..*. 

DO  50  1=1, K  -1 

IFC(N(I)  -  N(I+1)).NE.  0)  THEN 
FLAG=0 

ELSE 
END  IF 
50  CONTINUE 

PRINT  *,  'FLAG  IS: ' ,  FLAG 

***//  MAIN  PROGRAM  OUTER  LOOP  ST ART( EVERY  LOOP  IS  ONE  REPLICATION)//*** 


29 


ZFPREP  =  0. 

ZFAILS  =  0 
FAILS  =  0 
TOTREP  =  0 

_  _  "•*'  „  '1K*» 

LOOP  =0 

'10  IF(LOOP.LT.MAXREP)  THEN 
LOOP  =  LOOP  +  1 
IF(TOTREP.  LT.  MAXRUN)  THEN 
TOTREP  =  TOTREP  +  1 

SELCTA  =  1 
SELCTB  =  2 

***//  FILL  ARRAY  KEY(REPS)  WITH  INTEGERS  1  TO  K  TO  BE  USED  AS  OUTPUT 
***//  OF  THE  SUBROUTINE  SHSORT 
DO  95  REPS=1,  MAXREP 
KEYl(REPS)  =  REPS 
KEY2(REPS)  =  REPS 
KEY3(REPS)  =  REPS 
KEY4(REPS)  =  REPS 
95  CONTINUE 


***//  CALCULATE  NMAX  NOT  TO  PRINT  LONGER  THAN  THE  MAX  SAMPLE  SIZE 
***//  CALCULATE  THE  MAXIMUM  NO.  OF  TRIALS  AND  ITS  INDEX  NO.  //*** 


CALL  IMAX( N , K , NMAX , NINDX) 


***//  CALCULATE  THE  Ql’S  FROM  THE  GIVEN  MASTER  Q  AND  THE  Al’S 

'DO  115  1=1,  K 

QI( I)  =  MSTRQ  *  AI(I) 

115  CONTINUE 

DO  120  1=1,15 

DO  125  J=1 ,500 

UNIRV(I,J)  =999. 

TRIALSC I, J)  =  99999 
125  CONTINUE 

120  CONTINUE 

*** / /  DRAW  UNIFORM  (0,1)  RV'S  AND  CONVERT  TO  BERNOULLI  TRIALS  //*** 
DO  130  1=1,  K 

CALL  SRNDCSEED,  TEMP,  N(I),  MULT,  SORT) 

DO  135  J=l,  N( I) 

UNIRV( I , J)  =  TEMP(J) 

IF  (UNIRVC I , J). LE.  1  -  QIC  I) )  THEN 
TRIALSC I , J)  =  0 

ELSE 

TRIALSC I, J)  =  1 
END  IF 

135  CONTINUE 
130  CONTINUE 
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***//  CALCULATE  THE  NO."' OF  FAILURES  FOR  EACITCOMPONENT  //*** 

DO  150  1=1,  K 
FI( I)  =  0 
150  CONTINUE 

***//  CALCULATE  THE  F  SUB  I* S  AND  THE  GRAND  TOTAL  NO.  OF  FAILURES 
BIGF  =  0 
DO  155  1=1,  K 

DO  160  J=l,  N( I) 

FI(I)  =  FI(I)  +  TRIALS(I.J) 

160  CONTINUE 

IF(FICI).EQ.O)  THEN 

ZFPREP  =  ZFPREP  +  1 

ELSE 
END  IF 

***//  CALCULATE  THE  QHAT  SUB  I  *  S:  F  SUB  l'S  DIVIDED  BY  N  SUB  IS 
QHATI(I)  =  REAL(FKI))  /  N(I) 

BIGF  =  BIGF  +  FI( I) 

155  CONTINUE 

***//  CASE  WHERE  NO  COMPONENTS  HAVE  ANY  FAILURES  //*** 


IF(BIGF.EQ.O)  THEN 

ZFAILS  =  ZFAILS  +  1 
'""-DO  200  1=1,  K 

NIREAL(I)  =  REAL(N(I)) 

200  CONTINUE 

CALL  USMNMXC NIREAL ,K , INC , NIMIN ,NIMAX) 

DO  205  ALF=1,  MAXALF 

CALL  MDCHIC1  -  ALFA(ALF) ,2.  ,CHISQ(ALF, LOOP) .ERROR) 

RSHAT( ALF , LOOP)=  1  -  (CHISQ( ALF.LOOP)  /  REAL(2  *  NIMIN)) 
IFCFLAG.EQ.  1)  THEN 

RHTSTR( ALF, L00P)=1-(CHISQ( ALF, LOOP)  /  REAL( 2  *  N(l))) 

ELSE 
END  IF 

205  CONTINUE 

IF(PRNT. EQ.  1)  THEN 

"WRITE  ( 1,0007)  w'  “'**  ..  : 

WRITE(1,C1C15) 

WRITE( 1,3334)  QI 
WRITE( 1,0001) 

WRITEC 1 ,C1C15) 

DO  141  J=1 ,NMAX 

WRITE( 1,1111)  (UNIRVC I , J) ,  1=1,  K) 

141  CONTINUE 

WRITE( 1,0002) 

WRITEC 1,C1C15) 

DO  146  J=1 ,NMAX 

WRITEC 1,2222)  (TRIALSC I , J) ,  1=1,  K) 

146  CONTINUE 

WRITEC 1,0003) 

WRITEC 1,C 1C 15) 

WRITEC 1,3333)  FI 
WRITEC 1,0005) 

WRITEC 1,C 1C  15) 
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WRITE(1,3335)  N 
WRITEC 1,0004) 

WRITE( 1 ,C1C15) 

WRITEC 1,3334)  QHATI 

WRITE(  1,'C/'  'THE  MAXIMUM  Q  HAT  SUB  I  IS:”,  T40,  F8.5)’)  QHTMAX 
WRITE( 1, *  C/* 'THE  MAXI  Q  HAT  SUB  I  IS  ELMNT  NO. : f ' ,T40 , 15) ' )  QINDX 
WRITEC  1 , '  ( / "THE  GRAND  TOTAL  NO.  OF  FAILURES  IS:  ,T40,  15)')  BIGF 
ELSE 
END  IF 

DEGFRCLOOP)  =  2. 

GO  TO  10 

FT  CF 

FAILS  =  FAILS  +  1 
END  IF 

***//  FIND  THE  MAX  OF  THE  INDIVIDUAL  COMPONENT  UNRELIABILITIES 
CALL  RMAXC QHATI,  K,  QHTMAX,  QINDX) 

***//  PRINT  THE  RESULT  OF  THE  MAIN  OPERATING  ELEMENTS  OF  THE  PROGRAM 

IFCPRNT.EQ. 1)  THEN 
WRITEC 1,0007) 

WRITEC 1,C1C15) 

,  WRITEC  1,3334)  QI 
WRITEC 1,0001) 

WRITEC 1,C1C15) 

DO  140  J=1 ,NMAX 

WRITEC 1,1111)  CUNIRVCI , J) ,  1=1,  K) 

140  CONTINUE 

WRITEC 1,0002) 

WRITEC 1,C1C15) 

DO  145  J=1 ,NMAX 

WRITEC 1,2222)  (TRIALS(I, J) ,  1=1,  K) 

145  CONTINUE 

WRITEC 1,0003) 

WRITEC 1,C 1C 15) 

■“WRITEC  1,3333)  FI 
WRITEC 1,0005) 

WRITEC 1,C1C15) 

WRITEC 1,3335)  N 
WRITEC 1,0004) 

WRITEC 1 ,C1C15) 

WRITEC 1,3334)  QHATI 

WRITEC  1 ,'  C  /  THE  MAXIMUM  Q  HAT  SUB  I  IS:”,  T40,  F8.  5)')  QHTMAX 
WRITEC  1,'C/  "THE  MAXI  Q  HAT  SUB  I  IS  ELMNT  NO. :  f  ’  T40 , 15)  '  )  QINDX 
WRITEC  1,'C/ "THE  GRAND  TOTAL  NO.  OF  FAILURES  IS:  ,T40,  15)')  BIGF 
ELSE 
ENDIF 

***//  CALCULATE  THE  AHAT  SUB  I’S  C WEIGHT  ESTIMATES) 

SUMNAI  =  0. 

DO  165  1=1,  K 

AHATICI)  =  QHATIC)  /  QHTMAX 
SUMNAI  =  SUMNAI  +  N( I)  *  AHATI(I) 

165  CONTINUE 


32 


IF(PRNT.EQ. 1)  THEN 
WRITE( 1,0006) 

WRITE( 1 ,C1C15) 

WRITE( 1,3334)  AHATI 

ELSE 
END  IF 

***//  CALCULATE  1  REPLICATION  OF  UPPR  ALFA  C.  L.  ON  SYSTEM  RELIABILITY 
DEGFR(LOOP)  =  2  *  (1  +  BIGF) 

DO  170  ALF=1 ,  MAXALF 

CALL  MDCHK1  -  ALFA( ALF) ,DEGFR( LOOP) , CHISQC ALF, LOOP) ,  ERROR) 
QHTUPR( ALF, LOOP)  =  CHISQC ALF, LOOP)  /  (2  *  SUMNAI) 

IF(FLAG.EQ.  1)  THEN 

RHTSTR( ALF, LOOP)  =  1  -(CHISQC ALF, LOOP)  /  REAL(2*N(1))) 

ELSE 

:  v.  ,  END  IF 

*  +""  (ALF, LOOP),  ALFA  (ALF) 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  COMPNTS.  IN  SERIES 

CALL  RHTSRS(QHTUPR( ALF, LOOP),  AHATI, K,  RSHAT( ALF, LOOP ) ) 

*  +T40 ,F8.  5) ' )  RSHAT( ALF, LOOP) 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  BRIDGE  STRUCTURE  *** 

CALL  RHTBRG(QHTUPR(ALF, LOOP), AHATI ,K,RSHTBR( ALF, LOOP)) 

170  CONTINUE 

***//  THIS  ELSE  AND  ENDIF  ARE  FOR  THE  TEST  AGAINST  MAXRUN  *********** 

ELSE  *•**  —  -*■  "nw 

WRITEC 1 , ' (  "  ’  ’  /  ” PROGRAM  EXCEEDED  THE  MAX  NO.  OF  RUNS’ ' , 

+' '  ALLOWED  OF: ’ ' ,16) 1 )  TOTREP 
GOTO  9999 

END  IF 
GOTO  10 
END  IF 

***//  SORT  THE  ARRAYS  OF  SYSTEM  UNRELIABILITIES( 1  FOR  EACH  CONF.  LEVEL) 

DO  700  ALF=1 ,  MAXALF 

DO  800  REPS=1 ,  MAXREP 

TRANSQ(REPS)  =  QHTUPR( ALF, REPS) 

TRANSR(REPS)  =  RSHAT( ALF, REPS) 

TRNSTR(REPS)  =  RHTSTR( ALF, REPS) 

TRANBR(REPS)  =  RSHTBRC ALF, REPS) 

800  CONTINUE 

CALL  SHS0RT(TRANSQ,KEY1, MAXREP) 

CALL  SHS0RT(TRANSR,KEY2, MAXREP) 

CALL  SHS0RT(TRNSTR,KEY3 , MAXREP) 

CALL  SHS0RT(TRANBR,KEY4, MAXREP) 

DO  900  REPS=1 ,  MAXREP 

QHTUPRC ALF, REPS)  =  TRANSQ(REPS) 

RSHAT( ALF, REPS)  =  TRANSR(REPS) 

RHTSTRC ALF, REPS)  =  TRNSTR(REPS) 
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RSHTBRC  ALF , REPS )  =  TRANBR(REPS) 
900  CONTINUE 

700  CONTINUE 


***//  PRINT  OUTPUT  REPORT  HEADINGS  ************* 


WRITEC 1,6666) 

WRITEC 1,6667)  MAXREP 
WRITEC 1,6668)  K 
WRITE( 1,6669) 

IF(K.EQ.5)  THEN 

WRITEC 1,6699) 

ELSE 
END  IF 

WRITEC 1,6670)  MSTRQ 
WRITEC 1,6671) 

WRITEC 1.C1C15) 

WRITEC 1,3334)  AI 
WRITEC 1,0007) 

WRITEC 1.C1C15) 

WRITEC 1,3334)  QI 
WRITEC 1,0005) 

WRITEC 1.C1C15) 

WRITEC 1,3335)  N 
WRITEC 1,6674) 

WRITEC  2,6666)  —  — 

WRITEC 2, 6667)  MAXREP 
WRITEC 2, 6668)  K 
WRITEC 2,6669) 

IFCK. EQ. 5)  THEN 

WRITEC 1,6699) 

ELSE 
END  IF 

WRITEC 2,6670)  MSTRQ 
WRITEC 2,6671) 

WRITEC  2 ,C1C15) 

WRITEC 2, 3334)  AI 
WRITEC 2,0007) 

WRITEC  2 ,C1C15) 

WRITEC 2, 3334)  QI 
WRITEC 2,0005) 

WRITEC 2, C 1C 15) 

WRITEC 2, 3335)  N 
WRITEC 2,6674) 

WRITEC 2,' (’ ’SORTED  RSHAT  1  IS:  "  ,/10CF8. 5)) ' ) 
+CRSHATC 1 ,REPS) ,  REPS=1,  MAXREP) 

WRITEC  2,' ("SORTED  RSHAT  2  IS:  '  '  ,/10(F8.  5))  ) 
+(RSHAT(2,REPS) ,  REPS=1,  MAXREP) 

IFCFLAG.EQ.  1)  THEN 

WRITEC 2, '(' 'SORTED  RHTSTR  1  IS:  ,/10(F8. 5))  ) 
+( RHTSTR ( l.REPS) .  REPS=1,  MAXREP) 

WRITEC 2,' (f 'SORTED  RHTSTR  2  IS:  ,/10(F8. 5))  ) 
+( RHTSTR C 2 , REPS) ,  REPS=1,  MAXREP) 

ELSE 
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END  IF 

IF(K.EQ.  5)  THEN 

WRITE(2,'(  "SORTED  RSHTBR  1  IS:  "  ,/10(F8.  5))  1 ) 

+(RSHTBR( 1 , REPS) .  REPS=1,  MAXREP) 

WRITE(  2, '("SORTED  RSHTBR  2  IS:  "  ,/10(F8. 5)) ' ) 

+( RSHTBRC 2, REPS),  REPS=1,  MAXREP) 

ELSE 
END  IF 

***//  COMPUTE  THE  VALUE  RS  OF  THE  TRUE  SYSTEM  REL.  FNCTN. ( SERIES  SYSTEM) 
***//  AND  FOR  THE  5 -COMPONENT  BRIDGE  STRUCTURE 
CALL  RSRSCQI ,K.RS) 

WRITE(1,'("  ',////'  'THE  TRUE  SERIES  SYSTEM  ", 

+  ' 'RELIABILITY  VALUE  IS:  "  ,T51 ,F8. 5) ' )  RS 
CALL  RBRIDG(QI,K,RSBRDG) 

IF(K.  EQ.  5)  THEN 

WRITEC  1,‘C"  ",////"  THE  TRUE  BRIDGE  STRUCTURE  ", 

+  ' 'RELIABILITY  VALUE  IS:  "  ,T51,F8. 5) ' )  RSBRDG 
ELSE 
END  IF 

WRITEC 1,6675) 


***//  COMPUTE  THE  DIFFERENCE  'DELTAR'  BTWN^RS  AND  RSHAT  OF  THE  THEO 
****//  RETICAL  QUANTILE  GIVEN  BY  ALFA(MUST  USE  SORTED  RSHAT  ARRAY) 

IFCFLAG.EQ.  1)  THEN  ‘ 

WRITEC 1,5755) 

ELSE 
END  IF 

DO  450  ALF=1 ,  MAXALF 

QUANTL(ALF)  =  MAXREP  *  Cl  -  ALFACALF)) 

DELTARCALF)  =  RS  -  RSHATC ALF , QUANTLC ALF) ) 

DELBRGCALF)  =  RSBRDG  -  RSHTBRC ALF, QUANTLC ALF) ) 

IFCFLAG.EQ. 1)  THEN 

DELSTRC ALF)  =  RS  -  RHTSTR CALF, QUANTLC ALF)) 

WRITEC 1,5555)  MAXREP,  ALFACALF),  REALC  QUANTLC  ALF) ) 
WRITEC 1 , 5656 )  RHTSTR ( ALF .QUANTLC ALF) ) 

WRITEC 1,5657)  DELSTRC ALF) 

ELSE 
END  IF 

IFCK. EQ.  5)  THEN 

DELBRGCALF)  =  RSBRDG  -  RSHTBRC ALF, QUANTLC ALF) ) 

WRITEC 1,5555)  MAXREP,  ALFACALF),  REALC  QUANTLC  ALF) ) 
WRITEC 1,5666)  RSHTBRC ALF, QUANTLC ALF)) 

WRITEC 1,5667)  DELBRG(ALF) 

ELSE 
END  IF 

WRITEC 1,5555)  MAXREP,  ALFACALF),  REALC QUANTLC ALF ) ) 

WR ITE  C 1 , 5  5  5  6 )  RS H ATC  ALF , QUANTLC  ALF ) ) 

WRITEC 1,5557)  DELTAR (ALF) 

450  CONTINUE 

PRINT*,  'QUANTLC 1)  IS:',  QUANTLC 1) 

PRINT*,  'QUANTLC 2)  IS:',  QUANTLC 2) 

***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
j  j  ********  RSHAT  *******  j J-kidf 
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WRITE( 1,6676) 

DO  400  ALF=1 ,MAXALF 
TRUQNT(ALF)  =  0 

DO  500  REPS=1 ,  MAXREP 

-  -  u  DIFF(REPS)  =  RS  -  RSHAT( ALF, REPS) 

CONTINUE 

DO  600  REPS=1,  MAXREP 

IF( ABS(DIFF(REPS) ).  LE. EPS)  THEN 
TRUQNT(ALF)  =  REPS 


,F10. 5) * )  DIFF( 


WRITE( 1, '  (  ”  ' '  ,  / ' ' TRUE  CONFIDENCE  LIMIT  IS: 
F8.4)') 

(TRUQNT(ALF)  /  REAL( MAXREP) )  *  100. 

GO  TO  620 

ELSEIF(DIFF(REPS).LT.  0.  )  THEN 
TRUQNT(ALF)  =  REPS 
GO  TO  610 

ELSE 

*  END  IF  _ 

CONTINUE 

IF(TRUQNT(ALF).  EQ.  0.  )  THEN 
WRITE( 1,4443)  ALFA(ALF) 

WRITE(1,'(”  "./’’THE  SMALLEST”, 

”  DIFFERENCE  BETWEEN  RS  AND  RSHAT  IS:”,F10.5)r 
MAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1. )  THEN 
WRITE( 1,4442)  ALFA(ALF) 

WRITE(  1 , '  (  ”  ALL  RSHAT' '  , 

' '  ARE  GREATER  THAN  RS  ”  ) ’ ) 
ELSEIF(ABS(DIFF(TRUQNT(ALF))). LE. ABS(DIFF(TRUQNT( ALF) 
THEN 

WRITE ( 1,4444)  ALFA(ALF), 

L  (TRUQNT(ALF)  /  REAL( MAXREP))  *  100. 

WRITE ( 1 , 4445 )  RSHAT( ALF , TRUQNT( ALF) ) 

WRITEC 1,4446) 

ELSE 

WRITE( 1,4444)  ALFA(ALF) , 

l-  ((TRUQNT(ALF)-l)  /  RE AL( MAXREP) )  *  100. 

WRITE (1,4445)  RSHAT( ALF , TRUQNT( ALF ) - 1 ) 

WRITEC 1,4447) 

END  IF 
CONTINUE 


**//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
**j/  frtHWif****  RSHTBR  (BRIDGE)  *******  / /*** 


IF(K. EQ. 5)  THEN 
DO  401  ALF=1 .MAXALF 
TRUQNT(ALF)  =  0 

DO  501  REPS=1 ,  MAXREP 

DIFF(REPS)  =  RSBRDG  -  RSHTBR( ALF, REPS) 

CONTINUE 

DO  601  REPS=1,  MAXREP 

IF( ABS(DIFF(REPS) ). LE. EPS)  THEN 
TRUQNT(ALF)  =  REPS 

WRITEC  1, '(-!•• '  ”  ,/' 'TRUE* CONFIDENCE  LIMIT  IS:  ” 
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621  END 

401  CONTINUE 
ELSE 
END  IF 


F8.  4) * ) 

(TRUQNT(ALF)  /  REAL(MAXREP))  *  100. 

GO  TO  621 

ELSEIF(DIFF(REPS).  LT.  0.  )  THEN 
TRUQNTC ALF)  =  REPS 
GO  TO  611 

ELSE 
END  IF 
CONTINUE 

IF(TRUQNT( ALF).  EQ.  0.  )  THEN 
WRITEC 1,4443)  ALFA (ALF) 

WRITE(  1 , 1  (  "  ”,/ "THE  SMALLEST’  ’ , 

' '  DIFFERENCE  BETWEEN  RSBRDG  AND  RSHTBR  IS: ' ' , 
F10.  5 )••'-)  DIFF( MAXREP)  ^ 

ELSEIF(TRUQNTCALF).  EQ. 1.  )  THEN 
WRITEC 1,4442)  ALFA(ALF) 

WRITEC  1 , '  (  "  ",/'  'ALL  RSHTBR* '  , 

’ '  ARE  GREATER  THAN  RSBRDG' ' )' ) 

ELSEIF( ABS( DIFF( TRUQNT( ALF ) ) ) . LE. ABS(DIFF(TRUQNT( ALF) 
THEN 

WRITEC 1,4444)  ALFA(ALF) , 

(TRUQNT(ALF)  /  REAL(MAXREP) )  *  100. 

WRITEC 1 , 4449 )  RSHTBRC ALF .TRUQNTC ALF) ) 

WRITEC 1,4446) 

ELSE 

WRITEC 1,4444)  ALFA(ALF), 

C CTRUQNTC ALF) -1)  /  REALC MAXREP))  *  100. 

WRITEC  1,4449.)  RSHTBRC  ALF, TRUQNTC  ALF) -1) 

WRITEC 1,4447) 

END  IF 


***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
***  J  J  ********  RKTSTR  ********  //*** 


IFCFLAG.EQ.  1)  THEN 
DO  4400  ALF=1 , MAX ALF 
TRUQNTC ALF)  =  0 

DO  5500  REPS=1,  MAXREP 

DIFFCREPS)  =  RS  -  RHTSTRC ALF, REPS) 

5500  CONTINUE 

DO  6600  REPS=1,  MAXREP 

IFCABSCDIFFCREPS)). LE.EPS)  THEN 
TRUQNTC ALF)  =  REPS 
WRITEC 1 , ' C ' ’  ' ' , /  "  TRUE  CONFIDENCE  LIMIT  IS:  ' ' , 
+  F8.  4) ' ) 

+  CTRUQNTC ALF)  /  REALC MAXREP))  *  100 

GO  TO  6620 

ELSE IFCDIFFC REPS ).  LT.  0.  )  THEN 
TRUQNTC ALF)  =  REPS 
GO  TO  6610 

ELSE 
END  IF 
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CONTINUE 

IF(TRUQNT(ALF). EQ. 0. )  THEN 
WRITEC 1,4443)  ALFA(ALF) 

WRITEC 1 , ' (  ”  "  , /  "  THE  SMALLEST' ’ , 

' '  DIFFERENCE  BETWEEN  RS  AND  RHTSTR  IS:  "  , 

F9. 5) ' )  DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1.  )  THEN 
WRITE( 1,4442)  ALFA(ALF) 

WRITEC  1-r'  C  "  ’  '  ,*/ '  '  ALL  RHTSTR.'  '  , 

' '  ARE  GREATER  THAN  RS’ ' ) ’ ) 

ELSEIFCABSCDIFFCTRUQNT(ALF))).LE. ABSCDI1TCTRUQNTCALF)  -  1))) 
THEN 

WRITEC 1,4444)  ALFACALF), 

CTRUQNTCALF)  /  REALCMAXREP) )  *  100. 

WRITEC 1,4448)  RHTSTRC ALF , TRUQNTC ALF) ) 

WRITEC 1,4446) 

ELSE 

WRITEC 1,4444)  ALFACALF), 

C CTRUQNTCALF) -1)  /  REALCMAXREP ))  *  100. 

WRITEC 1,4448)  RHTSTRC  ALF , TRUQNTC  ALF )  - 1 ) 

WRITEC 1,4447) 

END  IF 
4400  CONTINUE 
ELSE 
END  IF 

***//  PRINT  THE  ARRAYS  PERTINENT  TO  THE  OUPUT  OF  EACH  REPLICATION  ***** 

IFCPRNT.EQ.  1)  THEN 
I  =  1 

185  WRITEC l.REPSHD)  ALFA(SELCTA) ,  ALFAC SELCTA) , 

+ALFACSELCTB) , ALFAC SELCTB) , ALFAC SELCTA) , ALFAC SELCTA) , ALFAC SELCTB) , 
+ALFACSELCTB) 

175  IFC-I.GE.  ( MAXREP  +  1))  THEN 
GOTO  180 

ELSE 

IFC  Cl*  EQ*  71).  OR.  Cl.  EQ.  211).  OR.  C  I.  EQ.  351).  OR.  C  I.  EQ.  491).  OR. 

+  CI-EQ.  631).  OR.  CI.EQ.  771).  OR.  Cl.  EQ.  911).0R.  CI.EQ.  1051)  )  THEN 
I  =  I  +  70 
WRITEC  1 , '  ( '  ’+")') 

GOTO  185 

ELSE 

WRITEC 1,3336)  I,  INTCDEGFR(I)) ,  CHISQ(1,I),  QHTUPR(1,I), 

+  CHISQ(2,I) ,  QHTUPRC 2,1) 

END  IF 

-IFC  I  +  70.  L&.  MAXREP y-’THEN  -•* 

WRITEC 1 ,3337)  1+70, INTCDEGFRC 1+70) ) .CHISQC 1,1+70) , 

+  QHTUPRC 1 , 1+70) , CHISQC  2 , 1+70) , QHTUPRC  2 , 1+70) 

ELSE 
END  IF 
1  =  1  +  1 
GOTO  175 
180  END  IF 
ELSE 
END  IF 


6600 
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9999  WRITE(  1 ,  *  (  *  *  THE  TOTAL  NO  OF  REPS  WAS:”,  18)’)  TOTREP 

WRITE ( 1 , ' ( ' 'THE  TOTAL  NO  OF  EFFECTIVE  REPS  WAS: ' '  18) ' )  LOOP 
WRITE (l.'C THE  TOTAL  NO  OF  NO  FAILURE  RUNS  WAS:'',  18)')  ZFAILS 
WRITE  (1, '(' 'AVERAGE  NO.  OF  COMPONENTS  PER  REPLICATION  WITH  ”, 

+'  'NO* FAILURES:  '-,F5.  2) '  )-ZFPREP  /  MAJ®EP 
WRITE( 1 , '  (  ' ' THE  TOTAL  NO  OF  RUNS  WITH  FAILURES  WAS:' ',18)')  FAILS' 
0008  FORMAT  (/  3X.’C  l',5X.'C  2', 

+5X,’C  3’,5X,fC  4',5X,fC  5’,5X,’C  6',5X,'C  7’ ,5X, 

+'C  8'  5X,'C  9 ' ,5X, ' C  10',4X,'C  11’ ,4X, 

+'C  12  ,4X, ' C  13',4X,'C  14'  4X,'C  15') 

0009  F0RMAT(/1X,'REP  NO' ,2X, ' DF  , IX. ' CHISQR( 1 ,F4.  3, ' ) ' , IX.  ^ 

+  ' QHTUPR( ’ ,F4.  3 , ' ) '  IX CHISQR( '  F4. 3 , ' ) '  IX, ' QHTUPRC ' ,F4.  3 , ' ) ' , 

+2X , ' REP  NO' , 2X , ' DF ' , IX  ' CHISQRC ' ,F4.  3 V , IX . 

+  ' QHTUPRC ' ,F4. 3 , ' ) ' , IX,  CHISQR( ' ,F4. 3,  ) ' . IX, ' QHTUPRC ’ ,F4.  3, ')' /) 
0001  FORMAT  C///'UNIFORM  RANDOM  DEVIATES  ARE: ') 

0002  FORMAT  C///' BERNOULLI  TRIALS  ARE:') 

0003  FORMAT  C///'TOTAL  NO.  OF  FAILURES  FOR  EACH  COMPONENT: ’) 

0004  FORMAT  (///'ESTIMATED  UNRELIABILITY  FOR  EACH  COMPONENT:') 

0005  FORMAT  (///'TOTAL  NUMBER  OF  MISSION  TESTS:  ') 

0006  FORMAT  C///' ESTIMATED  WEIGHTS  FOR  EACH  COMPONENT:') 

0007  FORMAT  C///'Q  I  FOR  EACH  COMPONENT: ') 

1111  FORMAT  C15F8.5) 

2222  FORMAT  (/IX, 15(I4,4X)) 

3333  FORMAT  C/1X, 15C I4.4X)) 

3334  FORMAT  (/15F8.5) 

3335  FORMAT  (/IX, 15 C 14, 4X)) 

3336  FORMAT  (T3 , 14, T9 , 13 ,T13 ,F11.  5 ,T27 ,F8. 5 ,T39 ,F11. 5 ,T53 ,F8. 5) 

3337  FORMAT  ( ’  +  ' ,T67 . 14 ,T73 , 13 ,T77 ,F11. 5 .T91 ,F8.  5 ,T103 ,F11.  5 ,T117 ,F8.  5) 

4442  FORMAT  ('  ',///'THE  RESULTING  (1  -  f ,F4. 3,’)  CONFIDENCE  , 

+’ LIMIT  IS: ' ,T50 , '  00.000  ') 

4443  FORMAT  ('  ',///'THE  RESULTING  (1  -  ' ,F4. 3,')  CONFIDENCE  , 

+  ' LIMIT  IS:  ' ,T50,' 100. 0000’ ) 

4444  FORMAT  ('  ',///' THE  RESULTING  (1  -  ’ ,F4. 3,')  CONFIDENCE  , 

+  'LIMIT  IS: ' ,T50,F8.  4) 

4445  FORMAT  ('  './'THE  RSHAT  VALUE  CLOSEST  TO  RS  IS:  '.T51.F8.  5)' 

4446  FORMAT  (’  './’(FIRST  NEGATIVE  DIFFERENCE)') 

4447  FORMAT  ('  './'(ELEMENT  PRECEEDING  FIRST  NEGATIVE  DIFFERENCE)  ) 

4448  FORMAT  ('  './'THE  RHTSTR  VALUE  CLOSEST  TO  RS  IS:  ’,T51,F8.5) 

4449  FORMAT  ('  './’THE  RSHTBR  VALUE  CLOSEST  TO  RSBRDG  IS:  ’.T51.F8.5) 

5555  FORMAT  (’  './//'THE  ' , 14, ' ( 1- ’ ,F4. 3 , ’ )  QUANTILE  IS:  '  .T49.F8.  3) 

5556  FORMAT  ('  './'THE  VALUE  OF  RSHAT  FOR  THAT  QUANTILE  IS:  ,T51,F8.5) 

5557  FORMAT  ('  './’THE  DIFFERENCES  -  RSHAT)  IS:  '  ,T51  ,F8.  5) 

5656  FORMAT  ('  ' ./'THE  VALUE  OF  RHTSTR  FOR  THAT  QUANTILE  IS: ’,T51,F8.S) 

5666  FORMAT  ('  './'THE  VALUE  OF  RSHTBR  FOR  THAT  QUANTILE  IS: ' ,T51,F8. 5) 

5657  FORMAT  ('  '  ,/'THE  DIFFERENCES  -  RHTSTR)  IS:  '  ,T51  ,F8.  5) 

5667  FORMAT  ('  './'THE  DIFFERENCES  -  RSHTBR)  IS:  '  ,T51,F8.  5) 

5755  FORMAT  ('  './//'SINCE  THE  NO.  OF  MISSION  TESTS  IS  THE  SAME  FOR  , 

+'  ALL  COMPONENTS  THE  CLOSED  FORM  SERIES  SYSTEM  RELIABILITY  ' , 

+ ”  '  RHTSTR ”  IS  COMPUTED '  )  ^ ;jl axxxaaj;xajuu.xaxxjiijuixjiajiax.i • 

+  ' Ir/rirlrtrlrirlr/rlrit  RUN  INPUT  SETTINGS  'MeMeMrkieMtlt-h'k 1  f 

4. 1  *  it -kirir -frir-kirir-ir-n iVVfVr iV Vr i'r i> Vr  '  ^ 

6667  FORMAT  ('  './/'NUMBER  OF  REPLICATIONS: ' ,T50,t4) 

6668  FORMAT  (’  './/'NUMBER  OF  COMPONENTS:  ' ,T50, 14) 

6669  FORMAT  ('  ’.//’SYSTEM  RELIABILITY  FUNCTION: ’ ,T50 ,' SERIES  ) 
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6699  FORMAT  ('  './/'SYSTEM  RELIABILITY  FUNCTION:  ' ,T50, ' BRIDGE' ) 

6670  FORMAT  ('  './/'MASTER  UNRELIABILITY  USED: ' .T50.F8. 5) 

6671  FORMAT  ('  './/'INPUT  VEIGHTS( A  SUB  I ' ' S) ; 

6674  FORMAT  ('+' , / / / ' *********************************************** ' j 
+ '**R  UN  RESULT  S*******************w********************' 

+  *  **** **********  1  ^ 

6675  FORMAT  ('  +  '  ,  / / /  '  ************************************************ ' 
+'  ESTIMATE  ERRORS  ********************************************* ' 

+ 1  ************** ' )  * 

6676  FORMAT  ('  +  ' . / // ' *********************************************** 1 
+'  TRUE  CONFIDENCE  LIMITS  *************************************** ^ 
*^v* ************ *  ^ 

END 


411 


APPENDIX  B.  INTERVAL  ESTIMATING  PROCEDURE  2 

PROGRAM  ZFNSCY 


TITLE:  BINOMIAL  INTERVAL  ESTIMATION  PROCEDURE  * 

ZERO  FAILURES  DISALLOWED;  WITH  SCALING  * 

AUTHOR:  E.  F.  BELLINI,  LT,  USN  * 

DATE:  NOV  89  * 

* 

THIS  PROGRAM  COMPUTES  THE  TRUE  CONFIDENCE  LEVEL  FOR  THE  ESTIMATE  * 
RELIABILITY  OF  A  SERIES  AND  BRIDGE  SYSTEM  GIVEN  THE  RELIABILITY  * 
OF  THEIR  COMPONENTS  * 

* 

IN  ITS  PRESENT  CONFIGURATION  THIS  PROGRlM  IS  SET  UP  TO  RUN  12  * 

TIMES  EACH  TIME  PRODUCING  1000  REPLICATIONS  USING  A  DIFFERENT  * 
SET  OF  INPUT  DATA.  RUN  THE  PROGRAM  FROM  CMS  BY  TYPING  'B1  EXEC'.* 
THE  REXX  EXEC  PROGRAM  * 

' B1 ’  CALLS  THE  INPUT  FILES  TO  BE  READ  AND  NAMES  THE  12  OUTPUT  * 
FILES  RESULTING  FROM  THE  12  CONSECUTIVE  RUNS.  BY  EDITING  THE  * 
INDEX  COUNTERS  I,  J,  K  OF  THE  ’ B1 ’  EXEC  ONE  CAN  RUN  ANY  USER-  * 
SPECIFIC  RUN  FROM  JUST  ONE  RUN  TO  ALL  12.  * 


VARIABLES  USED 


* 

AH  AT  I 

* 

AI 

* 

ALFA 

* 

BIGF 

Vc 

CHISQ 

V? 

COUNT 1 

V? 

C1C15 

V; 

DEGFR 

*.Y 

DELBRG 

* 

DELSTR 

Vr 

DELTAR 

Vr 

DIFF 

V: 

EPS 

yu 

ERROR 

Vf 

FAILS 

* 

FI 

jY 

FLAG 

V? 

INC 

* 

KEY1 

*>Y 

KEY2 

* 

KEY3 

* 

KEY4 

* 

KK 

* 

LOOP 

* 

MAXALF 

y? 

MAXREP 

Vr 

MAXRUN 

* 

MSTRQ 

WEIGHT  ESTIMATES  FOR  EACH  COMPONENT  * 

INPUT  WEIGHTS  FOR  EACH  COMPONENT  * 

LEVELS  OF  SIGNIFICANCE  * 

TOTAL  NO.  OF  FAILURES  FOR  EACH  REPLICATION  * 

CHI-SQUARE  RANDOM  VARIABLE  VALUE  * 

COUNTS  THE  NO.  OF  COMPONENTS  WITH  FAILURES 
FORMAT  LABEL  * 

DEGREES  OF  FREEDOM  * 

DIFFERENCE  FOR  BRIDGE  SYSTEM  * 

DIFFERENCE  FOR  SERIES  SYSTEM-  CLOSED  FORM  * 

DIFFERENCE  FOR  SERIES  SYSTEM  * 

DIFFERENCE  (TRUE  REL.  -  ESTIMATED  REL.  )  * 

SMALL  QUANTITYC CONSTANT)  * 

PARAMETER  FOR  IMSL  ROUTINE  * 

COUNTS  NO.  OF  REPLICATIONS  WITH  AT  LST.  1  FAILURE  * 
NO.  OF  FAILURES  FOR  EACH  COMPONENT ALL  MISSION  TST)* 
1  IF  ALL  COMP.  HAVE  SAME  NO.  OF  MISSION  TESTS  * 

INCREMENT  STEP  SIZE  FOR  ROUTINE  USMXMX  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  SIZING  PARAMETER  FOR  THE  MAX  NO  OF  COMPONENTS* 
COUNTS  NO.  OF  REPLICATION  PERFORMED  * 

MAX  NO.  OF  SIGNIFICANT  LEVELS  DESIREDC ARRAY  SIZING)* 
MAX  NO.  OF  REPLICATIONS  * 

MAX  NO.  OF  PROGRAM  ITERATIONS  ALLOWED  * 

MASTER  UNRELIAB1LITYCUSED  WITH  AI'S  TO  CALC.  QI'S)  * 
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k 

MULT 

MULTIPLIER  FOR  RANDOM  NO.  GENERATOR  SRND 

* 

k 

N 

NO.  OF  MISSION  TEST  FOR  EACH  COMPONENT  * 

IV 

SI  MAX 

MAX  NO.  OF  MISSION  TESTS 

k 

k 

NIMIM 

MIN  NO.  OF  MISSION  TESTS 

k 

* 

NINDX 

INDEX  NO.  OF  MAX  NO.  OF  MISSION  TESTS 

k 

* 

NTREAL 

NO.  OF  MISSION  TESTS  TRANSFORMED  TO  REAL 

k 

k 

NMAX 

MAX  NO.  OF  MISSION  TESTS  FOR  OUTPUT  CONTROL 

k 

k 

NPRNT 

FLAG  FOR  DETAILED  REPORT  OUTPUT 

k 

* 

PRNT 

SAME  AS  ABOVE (PARAMETER) 

k 

* 

QHATI 

UNRELIABILITY  ESTIMATES  FOR  EACH  COMPONENT 

k 

k 

QHTMAX 

LARGEST  QHATI 

k 

k 

QHTUPR 

UPPER  LIMIT  ON  SYSTEM  UNRELIABILITY 

k 

k 

QI 

INPUT  UNRELIABILIY  FOR  EACH  COMPONENT 

k 

k 

QINDX 

INDEX 

k 

k 

QUANTL 

QUANTILE 

k 

k 

REPSHD 

REPLICATIONS  HEADING  FORMAT  NUMBER 

k 

k 

RHTSTR 

SERIES  SYSTEM  RELIABILITY  ESTIMATE( CLOSED  FORM) 

k 

k 

FS 

TRUE  SERIES  SYSTEM  RELIABILITY 

k 

k 

RSBRDG 

TRUE  BRIDGE  SYSTEM  RELIABILITY 

k 

k 

RSHAT 

SERIES  SYSTEM  RELIABILITY  ESTIMATE 

k 

k 

RSHTBR 

BRIDGE  SYSTEM  RELIABILITY  ESTIMATE 

k 

k 

SEED 

PARAMETER 

k 

k 

SELCTA 

SIGNIFICANCE  LEVEL  SELECTION 

k 

k 

SELCTB 

SIGNIFICANCE  LEVEL  SELECTION 

k 

k 

SORT 

PARAMETER  FOR  ROUTINE  SRND 

k 

k 

SUMNAI 

SUM  OF  THE  PRODUCT  OF  Nl’S  AND  AI’S 

k 

i: 

TEMP 

TEMPORARY  ARRAY 

k 

k 

TOTREP 

TOTAL  NUMBER  OF  PROGRAM  ITERATIONS 

k 

k 

TRANBR 

TEMPORARY  ARRAY 

k 

k 

TRAN’SQ 

TEMPORARY  ARRAY 

k 

k 

TRAN SR 

TEMPORARY  ARRAY 

k 

k 

TRIALS 

BERNOULLI  TRIALS  ARRAY  (2-DIM) 

k 

k 

TRNSTR 

TEMPORARY  ARRAY 

k 

k 

TRUQNT 

TRUE  QUANTILE 

k 

k 

UNIRV 

UNIFORM  RANDOM  DEVIATES  (2-DIM) 

k 

k 

ZFAILS 

TOTAL  NUMBER  OF  REPLICATIONS  WITH  ZERO  FAILURES 

k 

•  r. 

k 

ZFPREP 

NO.  OF  COMPNTS .  WITH  ZERO  FAILURES  PER  REPLICATION 

k 

k 

Vf  *  *  Vr  *  yV  *  Vr  Vr  fr  fp  Vr  7V  Vr  *  Vr  fr  *  Vr  Vf  Vr  Vr  fr  tV  Vf  *  *  Yr  *  fr  Vf  TV  Vr  *  Vr  fc  *  TV  Vr  *  *  Vr  fr  A /r  tV  Vr  fr  tV  *  fr  fr  Vr  fr  Vr  fr  fr  Vr  Vr  Vr  Vr  tV  iV  *  * /r  * 


PARAMETER  (KK=10 ,MAXALF=2 ,NPRNT=0) 

PARAMETER  ( MAXREP=1000 ,  MAXRUN=10000,  EPS=. 000001) 

REAL*4  UNIRV( 15 ,1000) ,TEMP( 1000), QI(KK),  AI(KK) ,  AHATI(KK) 

REAL*4  QHATI(KK) ,  NMAX ,  QHTMAX ,  ALFA(MAXALF) 

RF.AL*4  DF(  5 ) , AALFA(  5 )  ,  5UMNAI ,  RSHAT(  MAXALF ,  MAXREP)  ,RS 
REAL*4  KEY  1  ( MAXREP) , KEY2( MAXREP) , KEY3 ( MAXREP ) ,TRNSTR( MAXREP) 
REAL*4  DEGFRC MAXREP),  QHTUPR( MAXALF, MAXREP) , CHI SQ( MAXALF , MAXREP) 
REAL* 4  RHTSTR( MAXALF , MAXREP) 

REAL* 4  DELTARC MAXALF) ,  TRANSQ( MAXREP) , TRAN SR (MAXREP) .DIFF(MAXREP) 
REAL*4  DELSTRC MAXALF) ,  NN(KK) 

REAL* 4  RSHTBR( MAXALF, MAXREP) ,DELBRG( MAXALF) ,KEY4( MAXREP) 

REAL*4  TRANBRC MAXREP ) ,  RSBRDG 

REAL* 4  SCALES,  MINQI ,  MAXQI ,  ZFPREP.MSTRQ 


INTEGER  SEED,  MULT,  SORT,  TRIALS( 15 , 1000) ,  BIGF ,  FI(KK) ,  N(KK) 


r 


INTEGER  NINDX,  QINDX,  ERROR,  REPS,  SELCTA,  SELCTQ ,  TOTREP 
INTEGER  COUNT1,  C1C15,  REPSHD,  SELCTB ,  ALF,  FLAG,  LOOP,PRNT 
INTEGER  QUANTL(MAXALF) ,  TRUQNT(MAXALF) ,  ZFAILS,  FAILS 
INTEGER  INC 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


DATA  SEED/123457/,  MULT/1/,  INC/1/ 

DATA  AALFA/ .  01 , .  05 ,.  9 , .  95 , .  99/ ,  DF/1,5,10,30,40/ 

DATA  ALFA/. 20,.  050/ 

DATA  MSTRQ/O.  1,  .05/ 

DATA  AI/.  20,.  10,.  05,.  05,.  05/ 

DATA  AI/.  20,.  10,.  05,.  05,.  05,.  05,.  05,.  05,.  05,.  05/ 

DATA  AI/.  05,.  05,.  05,.  05,.  05,.  05,.  05,.  05,.  05,.  05/ 

DATA  AI/.  05 , .  05 , .  05 , .  05 , .  05 , .  05  , .  05 , .  05  , .  05  , .  05  , .  05 , .  05 , .  05 , 

+.05,.  05/ 

DATA  AI/.  1,.  1,.  1,.  1,1,.  05,.  05,.  2,.  25,.  5/ 

DATA  AI/.  1,.  1,1.  ,.05,. 2/ 

DATA  N/100,100,100,100,100/.. 

DATA “N/10, 6, 5, 2, 10, 4, 5, 6, 7, 2/ 

DATA  N/50, 50, 50, 50, 50, 50, 50, 50, 50, 50/ 

DATA  N/ 100 , 100 , 100 , 100 ,100,100, 100 ,100,100, 100/ 

DATA  N/300, 300, 300, 300, 300, 300, 300, 300, 300, 300/ 

DATA  N/20, 20, 20, 20, 20, 20, 20, 20, 20, 20/ 

DATA  N/ 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 
+1000,1000,1000,1000/ 

DATA  SORT/O/ 


ASSIGN  8  TO  C1C15 
ASSIGN  9  TO  REPSHD 

CALL  EXCMS( 'FILEDEF  01  DISK  B  OUTPUT  AI  (LRECL  132  ') 
CALL  EXCMSC FILEDEF  02  DISK  JNK  OUTPUT  AI  (LRECL  132  ') 

SELCTQ  =  1 
PRNT  =  NPRNT 


DO  12  1=1, KK 

AI(I)  =  9999. 

NCI  5  =  99999999 
12  CONTINUE 

RE  AD  ( 0  3 ,  *  )  K ,  MSTRQ 

DO  11  1=1, K 
READC03,*)  AI( I) ,N( I) 
11  CONTINUE 


***//  INITIALIZE  THE  QHTUPR  ARRAY  OF  UNRELIABILITY  REPLICATIONS,  //*** 

*  RSHAT  ARRAY  OF  ESTIMATE  SYSTEM  RELIABILITY  REPLICATIONS  * 

*  AND  RHTSTR  ARRAY  OF  EST.  SYST.  REL.  FOR  A  SERIES  SYST  WHEN  * 

***//  ALL  THE  COMPONENT  MISSION  TESTS  ARE  EQUAL  IN  NUMBER  //*** 


DO  172  ALF=1,MAXALF 

DO  173  REPS=1 ,MAXREP 

QHTUPRC ALF, REPS)  =  0. 
RSHAT( ALF, REPS)  =  0. 
RHTSTRC ALF, REPS)  =  0. 


I 

I 


RSHTBR(ALF,REPS)  =  0. 

173  CONTINUE  —  ■*- 

172  CONTINUE  ' 

***//  SET  FLAG  TO  1  IF  ALL  COMPONENTS  HAVE  SAME  NO.  OF  MISSION  TESTS**** 
FLAG=1 

DO  50  1=1 ,K  -1 

IF( (N( I)  -  N( 1+1) ) .  NE.  0)  THEN 
FLAG=0 

ELSE 
END  IF 
50  CONTINUE 

PRINT  *,  ’FLAG  IS: 1  ,  FLAG 

***//  MAIN  PROGRAM  OUTER  LOOP  START( EVERY  LOOP  IS  ONE  REPLICATION)//*** 

ZFPREP  =  0. 

FAILS  =  0 
ZFAILS  =  0 
COUNT 1  =  0 
TOTREP  =  0 
LOOP  =  0 

10  IF(LOOP.  LT.MAXREP)  THEN 

IFCTOTREP.  LT.MAXRUN)  THEN 
SELCTA  =  1 
SELCTB  =  2 

***//  FILL  ARRAY  KEY(REPS)  WITH  INTEGERS  1  TO  K  TO  BE  USED  AS  OUTPUT 
***//  OF  THE  SUBROUTINE  SHSORT 

DO  95  REPS=1,  MAXREP 
KEYl(REPS)  =  REPS 
KEY2(REPS)  =  REPS 
KEY3(REPS)  =  REPS 
KEY4CREPS)  =  REPS 
95  CONTINUE 

***//  CALCULATE  NMAX  NOT  TO  PRINT  LONGER  THAN  THE  MAX  SAMPLE  SIZE  //*** 
***//  CALCULATE  THE  MAXIMUM  NO.  OF  TRIALS  AND  ITS  INDEX  NOV  //*** 

CALL  IMAX( N ,  K , NMAX , NINDX) 


***//  CALCULATE  THE  Ql’S  FROM  THE  GIVEN  MASTER  Q  AND  THE  Al’S 

DO  115  1=1,  K 

QIC  ID  =  MSTRQ  *  AI(I) 

115  CONTINUE 


DO  120  1=1,15 

DO  125  J=1 ,500 

UNIRVC I , J)  =  999. 
TRIALS(I.J)  =  99999 
125  CONTINUE 
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120  CONTINUE 

***//  CALCULATE  THE  SCALING  FACTOR  'SCALEN' 

CALL  USMNMX( QI , K , INC , MINQI , MAXQI ) 

SCALEN  =  MAXQI  /  MINQI 

***//  draw  UNIFORM  (0,1)  RV’S  AND  CONVERT  TO  BERNOULLI  TRIALS 
DO  130  1=1,  K 

CALL  SRND(SEED,  TEMP,  N(I),  MULT,  SORT) 

DO  135  J=l,  N( I) 

UNIRV( I , J)  =  TEMP(J) 

IF  (UNIRV( I , J).  LE.  1  -  QIC  I) )  THEN 
TRIALS( I , J)  =  0 

ELSE 

TRIALS( I , J)  =  1 
END  IF 

135  CONTINUE 
130  CONTINUE 


***//  CALCULATE  THE  NO.  OF  FAILURES  FOR  EACH  COMPONENT 

DO  150  1=1,  K 
FI( I )  =  0 
150  CONTINUE 

***//  CALCULATE  THE  F  SUB  I'S  AND  THE  GRAND  TOTAL  NO.  OF  FAILURES 

BIGF  =  0 
DO  155  1=1,  K 

DO  160  J=l,  N( I) 

FI( I)  =  FI( I)  +  TRIALS( I , J) 

160  —  CONTINUE*"  — 1  -*» 

***//  CALCULATE  THE  QHAT  SUB  I'S:  F  SUB  I'S  DIVIDED  BY  N  SUB  I'S 

IF(FI( I) .  EQ.  0)  THEN 

QHATI(I)  =  1.  /  (SCALEN  *  N(I)) 

ZFPREP  =  ZFPREP  +  1 

ELSE 

QHATI(I)  =  REAL(FIU))  /  N(I) 

END  IF 

BIGF  =  BIGF  +  FI( I ) 

155  CONTINUE 

***//  FIND  THE  MAX  OF  THE  INDIVIDUAL  COMPONENT  UNRELIABILITIES 

CALL  RMAX(QHATI ,  K,  QHTMAX ,  QINDX) 

IF(PRNT.EQ.  1)  THEN 
WRITE( 1,0007) 

WRITE( 1 ,C1C15) 

VRITEC 1,3334)  QI 
WRITE( 1,0001) 


j  /***<■ ••>»- 


/ /*** 


/ /Irk* 


f  /*** 


//*** 
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WRITE(1,C1C15) 

DO  140  J=1 ,NMAX 

WRITE( 1,1111)  (UNIRV( I » J) >  1=1,  K) 

140  CONTINUE 

WRITE( 1,0002) 

WRITE( 1 ,C1C15) 

DO  145  J=1 ,NMAX 

WRITE( 1,2222)  (TRIALSC I , J) ,  1=1,  K) 

145  CONTINUE 

WRITE( 1,0003) 

WRITE( 1 ,C1C15) 

WRITE( 1,3333)  FI 
WRITE( 1,0005) 

WRITE( 1 ,C1C15) 

WRITE( 1,3335)  N 
WRITE( 1,0004) 

WRITE( 1 ,C1C15) 

,,WRITE(  1.3334)  QHATI 

WRITE  ( 1 , '  (  / '  '  THE  MAXIMUM  Q  HAT  SUB  I  IS:”,  T40.  F8.5)')  QHTMAX 
WRITE(1,'(/' 'THE  MAXI  Q  HAT  SUB  I  IS  ELMNT  NO. : 1 '  T40 , 15) ' )  QINDX 
WRITE(1,'(/' ’THE  GRAND  TOTAL  NO.  OF  FAILURES  IS:  ,T40,  15)')  BIGF 
ELSE 
ENDIF 

***//  TEST  FOR  A  REP  WITH  AT  LST  TWO  COMP  WITH  AT  LST  ONE  FAILURE  EACH 


162 


DO  162  J=l,  K 
NN( J)  =  0 

IF(FI( J). GE. 1)  THEN 

COUNT 1  =  COUNT 1  +  1 
NN(  J)=  NXJ) 

ELSE 

NN( J)  =  N(J)  *  SCALEN 
END  IF 
CONTINUE 

IF ( COUNT l.EQ.  0)  THEN 

ZFAILS  =  ZFAILS  +  1 


ELSE 

FAILS  =  FAILS  +  1 
END  IF 

IFCCOUNTl.GE. 2)  THEN 
LOOP  =  LOOP  +  1 
TOTREP  =  TOTREP  +  1 

ELSE 

TOTREP  =  TOTREP  +  1 
GO  TO  10 
END  IF 


***//  CALCULATE  THE  AHAT  SUB  I'S  (WEIGHT  ESTIMATES)  //*** 

SUMNAI  =  0. 

DO  165  1=1,  K 

AHATI(I)  =  QHATI(I)  /  QHTMAX 

*  +(FI(L),L=1,K) 

*  +(NN(L),L=1,K) 

SUMNAI  =  SUMNAI  +  NN(I)  *  AHATI(I) 
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165  CONTINUE 

IF(PRNT.  EQ.  1)  THEN 
WRITE( 1,0006) 

WRITE( 1 ,C1C15) 

WRITEC 1,3334)  AHATI 

ELSE 
END  IF 

***//  CALCULATE  1  REPLICATION  OF  UPPR  ALFA  C.  L.  ON  SYSTEM  RELIABILITY 
DEGFR(LOOP)  =  2  *  (1  +  BIGF) 

DO  170  ALF=1 ,  MAXALF 

CALL  MDCHIU  -  ALFA( ALF) ,DEGFR(LOOP) ,CHISQ( ALF, LOOP) ,  ERROR) 
QHTUPRC ALF , LOOP)  =  CHISQ(ALF.LOOP)  /  (2  *  SUMNAI) 

IF(FLAG.EQ.  1)  THEN 

RHTSTRC ALF, LOOP)  =  1  -( CHI SQ( ALF, LOOP)  /  REAL(2*N(1))) 

ELSE 
END  IF 

*  +  (ALF, LOOP),  ALFA(ALF) 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  COMPNTS.  IN  SERIES 

CALL  RHTSRS(QHTUPR( ALF, LOOP),  AHATI, K,  RSHAT( ALF, LOOP) ) 

*  +T40 ,F8.  5) ' )  RSHAT( ALF, LOOP) 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  BRIDGE  STRUCTURE  *** 
CALL  RHTBRGC QHTUPRC ALF, LOOP), AHATI, K,RSHTBR( ALF, LOOP)) 

170  CONTINUE 

***//  THIS  ELSE  AND  END IF  ARE  FOR  THE  TEST  AGAINST  MAXRUN  //*** 

ELSE 

WRITE(  !,*(’’  '  '  /  ”  PROGRAM  EXCEEDED  THE  MAX  NO.  OF  RUNS’  '  , 

+  ’ '  ALLOWED  OF:  ' '  ,16) f)  TOTREP 
GOTO  9999 

END  IF 
GOTO  10 
END  IF 

***//  SORT  THE  ARRAYS  OF  SYSTEM  UNRELIABILITIESC 1  FOR  EACH  CONF.  LEVEL) 

DO  700  ALF=1 ,  MAXALF 

DO  800  REPS=1 ,  MAXREP 

TRANSQCREPS)  =  QHTUPRC ALF, REPS) 

TRANSR(REPS)  =  RSHAT( ALF, REPS) 

TRNSTR(REPS)  =  RHTSTRC ALF, REPS) 

TRANBR(REPS)  =  RSHTBRC ALF, REPS) 

800  CONTINUE 

CALL  SHSORTCTRANSQ, KEY 1, MAXREP) 

CALL  SHSORTC TRANSR , KEY2 , MAXREP ) 

CALL  SHSORTC TRNSTR , KEY3 , MAXREP ) 

CALL  SHSORTC TRANBR ,KEY4 , MAXREP) 

DO  900  REPS=1 ,  MAXREP 
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QHTUPR( ALF , REPS )  =  TRANSQ(REPS) 
RSHAT(ALF,REPS)  =  TRANSR(REPS) 
RHTSTR( ALF , REPS )  =  TRNSTR(REPS) 
RSHTBRC ALF, REPS)  =  TRANBR(REPS) 
900  ...  CONTINUE. 

700  CONTINUE 


***//  PRINT  OUTPUT  REPORT  HEADINGS 


WRITEC 1,6666) 

WRITEC 1,6667)  MAXREP 
WRITE( 1,6668)  K 
WRITEC 1,6770)  SCALEN 
WRITEC 1,6669) 

IFCK.EQ.  5)  THEN 

WRITEC 1,6699) 

ELSE 
END  IF 

WRITEC 1,6670)  MSTRQ 
WRITEC 1,6671) 

WRITEC 1,C1C15) 

WRITEC 1,3334)  AI 
WRITEC 1,0007) 

WRITEC 1,C1C15) 

WRITEC 1,3334)  QI 
WRITEC 1,0005) 

WRITEC 1.C1C15) 

WRITEC 1,3335)  N 
...WRITEC  1,6674) 

WRITEC 2, 6666) 

WRITEC 2, 6667)  MAXREP 
WRITEC 2,6668)  K 
WRITEC 2,6770)  SCALEN 
WRITEC 2, 6669) 

IFCK.EQ.  5)  THEN 

WRITEC 1,6699) 

ELSE 
END  IF 

WRITEC 2,6670)  MSTRQ 
WRITEC 2, 6671) 

WRITEC 2, C1C15) 

WRITEC 2, 3334).  AI 
WRITEC 2, 0007) 

WRITEC 2, C 1C 15) 

WRITEC 2, 3334)  QI 
WRITEC 2, 0005) 

WRITE(2,C1C15) 

WRITEC 2, 3335)  N 
WRITEC 2, 6674) 

WRITEC2,'C  "SORTED  RSHAT  1  IS:" 
+CRSHATC 1,REPS) ,  REPS=1,  MAXREP) 
WRITEC  2 , '  C  "  SORTED  RSHAT  2  IS:  " 
+C RSHATC  2 , REPS ) ,  REPS=1,  MAXREP) 
IFCFLAG.EQ.  1)  THEN 

WRITEC  2 , ' (  "  SORTED  RHTSTR  1 


,/15CF8.5))') 

,/15CF8.5))') 


IS:  ' '  ,/10(F8. 5)) ' ) 


- 


I /*** 


■ 


4S 


+(RHTSTR(1,REPS).  REPS=1,  MAXREP) 
WRITE(2,'(  "SORTED  RHTSTR  2 
+(RHTSTR( 2 , REPS) ,  REPS=1,  MAXREP) 
ELSE 
END  IF 

IFCK.EQ.5)  THEN 

VRITE(2,'(  "SORTED  RSHTBR  1 
+(RSHTBR( 1 , REPS) .  REPS=1,  MAXREP) 
WRITE(  2,' ("SORTED  RSHTBR  2 
+(RSHTBR(2,REPS) ,  REPS=1,  MAXREP) 
ELSE 
END  IF 


IS:  ' ' ,/10(F8. 5)) ' ) 

IS:  '  '  ,/10(F8.5))') 
IS:  ' ' ,/10(F8. 5) ) 1 ) 


***//  COMPUTE  THE  VALUE  RS  OF  THE  TRUE  SYSTEM  REL.  FNCTN. (SERIES  SYSTEM) 
***//  AND  FOR  THE  5 -COMPONENT  BRIDGE  STRUCTURE 


CALL  RSRS(QI.K.RS) 

WRITE(1,'("  ",////'  'THE  TRUE  SERIES  SYSTEM  ", 

+  ' ' RELIABILITY  VALUE  IS:  "  ,T51,F8.  5) ' )  RS 
CALL  RBRIDG(QI,K,RSBRDG) 

IF(K. EQ.  5)  THEN 

WRITE(1,’("  ".////"THE  TRUE  BRIDGE  STRUCTURE  ", 
+  "  RELIABILITY  VALUE  IS:  "  ,T5 1  ,F8.  5 ) '  )  RSBRDG 
ELSE 
END  IF 

WRITE( 1,6675) 


***//  COMPUTE  T51E  DIFFERENCE  'DELTAR'  BTWN.  RS  AND  RSHAT  OF  THE  THEO 
***//  RETICAL  QUANTILE  GIVEN  BY  ALFA(MUST  USE  SORTED  RSHAT  ARRAY) 

IF(FLAG.  EQ.  1)  THEN 
WRITE( 1,5755) 

ELSE 
END  IF 

DO  45 0  ALF=1 ,  MAXALF 

QUANTL(ALF)  =  MAXREP  *  (1  -  ALFA(ALF)) 

DELTAR(ALF)  =  RS  -  RSHAT( ALF , QUANTL( ALF) ) 

DELBRG(ALF)  =  RSBRDG  -  RSHTBR( ALF , QUANTL( ALF ) ) 

IF(FLAG.  EQ.  1)  THEN 

DELSTR(ALF)  =  RS  -  RHTSTR ( ALF , QUANTL( ALF ) )  . 

WRITE( 1,5555)  MAXREP,  ALFA(ALF),  RE AL( QUANTL( ALF ) ) 
WRITE(  1 , 5  65  6 )  RHTSTR ( ALF , QUANTL( ALF) ) 

WRITE( 1,5657)  DELSTR(ALF) 

ELSE 
END  IF 

IF(K. EQ.  5)  THEN 

DELBRG(ALF)  =  RSBRDG  -  RSHTBR( ALF , QUANTL( ALF) ) 

VRITE( 1,5555)  MAXREP,  ALFA(ALF),  REAL( QUANTL( ALF) ) 

WR ITE ( 1 , 5  6  6  6 )  RSHTB  R ( ALF , QUANTL( ALF ) ) 

WRITE( 1,5667)  DELBRG(ALF) 

ELSE 
END  IF 

VRITE( 1,5555)  MAXREP,  ALFA( ALF) ? REAL(QUANTL( ALF) ) 

WRITE( 1,5556)  RSHAT( ALF ,QUANTL( ALF) ) 

VRITE( 1,5557)  DELTAR(ALF) 

450  CONTINUE 
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PRINT*,  ’QUANTL(l)  IS:',  QUANTL(l) 
PRINT*,  ' QUANTLC 2 )  IS:',  QUANTL(2) 


***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
*** / /  ********  RSHAT  *******  j /*** 


WRITEC 1,6676)  '  w 

DO  400  ALF=1 ,MAXALF 
TRUQNT(ALF)  =  0 

DO  500  REPS=1,  MAXREP 

DIFF(REPS)  =  RS  -  RSHAT(ALF.REPS) 

500  CONTINUE 

DO  600  REPS=1,  MAXREP 

IF(ABS(DIFF(REPS) ).  LE.  EPS)  THEN 
TRUQNTC ALF)  =  REPS 

WRITE( 1 , ' (  "  "  ,  /  "  TRUE  CONFIDENCE  LIMIT  IS: ' ' , 

+  F8.  4)  '  ) 

+  (TRUQNT(ALF)  /  RE AL( MAXREP))  *  100. 

GO  TO  620 

ELSEIF(DIFF(REPS).  LT.  0.  )  THEN 
TRUQNT(ALF)  =  REPS 
GO  TO  610 

ELSE 
END  IF 

600  CONTINUE 

610  1F(TRUQNT(ALF).  EQ. 0. )  THEN 

WRITE ( 1,4443)  ALFA(ALF) 

WRITE(  1 , '  (  "  * './"THE  SMALLEST", 

+  "  DIFFERENCE  BETWEEN  RS  AND  RSHAT  IS: ' '  ,F10. 5) ' )  DIFF( 

+  MAXREP) 

ELSEIFCTRUQNT(ALF).  EQ.  1.  )  THEN 
WRITEC 1,4442)  ALFA(ALF) 

WRITE(1,'("  ALL  RSHAT", 

+  "  ARE  GREATER  THAN  RS")’) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.  ABS(DIFF(TRUQNT( ALF)  -  1))) 

+  THEN 

WRITEC 1,4444)  ALFA (ALF) , 

+  (TRUQNT(ALF)  /  REAL( MAXREP) )  *  100. 

WRITEC 1,4445)  RSHAT ( ALF , TRUQNTC ALF) ) 

WRITEC 1,4446) 

ELSE 

WRITEC 1,4444)  ALFA(ALF), 

+  (( TRUQNTC ALF) -1)  /  REALC MAXREP))  *  100. 

WRITEC 1,4445)  RSHAT( ALF , TRUQNTC ALF) -1) 

WRITEC 1,4447) 

620  END  IF 

400  CONTINUE 


***//  find  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
*** j j  ********  RSHTBR  (BRIDGE)  *******  / /*** 

IFCK.EQ. 5)  THEN 
DO  401  ALF=1 .MAXALF 
TRUQNTC ALF)  =  0 

DO  501  REPS=1 ,  MAXREP 
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DIFF(REPS)  =  RSBRDG  -  RSHTBRC ALF, REPS) 

501  CONTINUE 

DO  601  REPS=1 ,  MAXREP 

IF(ABS(DIFF(REPS)).LE.  EPS)  THEN 
TRUQNT(ALF)  =  REPS 

WRITEC 1 , ' ( ' '  "  , /  "  TRUE  CONFIDENCE  LIMIT  IS:  ”  , 

+  F8.  4)  '  ) 

+  (TRUQNT(ALF)  /  REAL( MAXREP) )  *  100. 

GO  TO  621 

ELSEIF(DIFF(REPS). LT. 0. )  THEN 
TRUQNT(ALF)  =  REPS 
GO  TO  611 

ELSE 
END  IF 

601  CONTINUE 

611  IF(TRUQNT(ALF).  EQ.  0.  )  THEN 

WRITE( 1,4443)  ALFA(ALF) 

WRITEC 1 , ' C  "  "  , /  "  THE  SMALLEST' ' , 

+  ,  .  '  '  DIFFERENCE  BETWEEN  RSBRDG  AND  RSHTBR  IS:  ”  , 

*  F10.5)')  DIFFCMAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1.  )  THEN 
WRITEC 1,4442)  ALFA(ALF) 

WRITEC  1,'("  ",/’  'ALL  RSHTBR"  , 

+  "ARE  GREATER  THAN  RSBRDG' ' )' ) 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.  ABSCDIFFCTRUQNTCALF)  -  1))) 
+  THEN 

WRITEC 1,4444)  ALFACALF) > 

+  CTRUQNTCALF)  /  REALC MAXREP))  *  100. 

WRITEC 1 , 4449 )  RSHTBRC ALF ,TRUQNTC ALF) ) 

WRITEC 1,4446) 

ELSE 

.  WRITEC  U4444)  ALFACALF),  ^ 

+  C C TRUQNTC ALF) -1)  /  REALC MAXREP))  *  100. 

WRITEC 1,4449)  RSHTBR C ALF, TRUQNTC ALF) -1) 

WRITEC 1,4447) 

621  END  IF 

401  CONTINUE 
ELSE 
END  IF 


***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
***  J I  •irMeMeMrtt  RHTSTR  ********  J  /*** 


IFCFLAG.EQ.  1)  THEN 
DO  4400  ALF=1 .MAXALF 
TRUQNTC ALF)  =  0 

DO  5500  REPS=1 ,  MAXREP 

DIFFCREPS)  =  RS  -  RHTSTRC ALF, REPS) 

5500  CONTINUE 

DO  6600  REPS=1 ,  MAXREP 

IFCABS(DIFFCREPS)).  LE.  EPS)  THEN 
TRUQNTC ALF)  =  REPS 

WRITEC  1,’C"  ’ '  ,  /  '  '  TRUE  CONFIDENCE  LIMIT  IS:" 
+  F8.4V) 

+  CTRUQNTCALF)  /  REALC MAXREP) )  *  lOO. 
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GO  TO  6620 

ELSEIF(DIFF(REPS). LT.  0. )  THEN 
TRUQNT(ALF)  =  REPS 
GO  TO  6610 

ELSE 
END  IF 

6600  CONTINUE 

6610  IF(TRUQNT(ALF).  EQ.  0. )  THEN 

WRITE( 1,4443)  ALFA(ALF) 

WRITEC 1 , ’ (  ”  "  , /  "  THE  SMALLEST’ ' , 

+  ' '  DIFFERENCE  BETWEEN  RS  AND  RHTSTR  IS:  "  , 

+  F9. 5) ' )  DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1.  )  THEN 
WRITEC 1,4442)  ALFA(ALF) 

WRITEC  1,’("  ALL  RHTSTR' ’  , 

+  "ARE  GREATER  THAN  RS  "  ) ' ) 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.ABSCDIFF(TRUQNT(ALF)  -  1))) 

+  THEN 

WRITEC 1,4444)  ALFACALF), 

+  CTRUQNTCALF)  /  REALC MAXREP))  *  100. 

WRITEC 1,4448)  RHTSTRC ALF , TRUQNTC ALF) ) 

WRITEC 1,4446) 

ELSE 

WRITEC 1,4444)  ALFACALF), 

+  C (TRUQNTC ALF) -1)  /  REALC MAXREP) )  *  100. 

WRITE (1,4448)  RHTSTRC  ALF .TRUQNTC  ALF) - 1 ) 

WRITEC 1,4447) 

6620  END  IF 

4400  CONTINUE 
ELSE 
END  IF 

***//  print  THE  ARRAYS  PERTINENT  TO  THE  OUPUT  OF  EACH  REPLICATION  //*** 

IF(PRNT. EQ. 1)  THEN 
I  =  1 

185  WRITEC l.REPSHD)  ALFA(SELCTA) ,  ALFA( SELCTA) , 

+ALFA(SELCTB) ,ALFA( SELCTB) ,ALFA( SELCTA) ,ALFA( SELCTA) ,ALFA(SELCTB) , 
+ALFACSELCTB) 

175  IF( I.  GE. (MAXREP  +  1))  THEN 
■GOjlO  180 

ELSE 

IF(  (I.EQ.  71).  OR.  (  I.  EQ.  211) .  OR.  ( I.  EQ.  351).  OR.  ( I.  EQ.  491).  OR. 

+  (I.  EQ.  631).  OR.  (  I.  EQ.  771).  OR.  (I.  EQ.  911).  OR.  (I.  EQ.  1051)  )  THEN 

I  =  I  +  70 
WRITEC  1, '("  +  ")’) 

GOTO  185 

ELSE 


WRITEC 1,3336)  I,  INT(DEGFR( I)) ,  CHISQ(1,I),  QHTUPR(l.I), 
+  CHISQC2.I),  QHTUPRC 2,1) 

END  IF 

IF( I  +  70.  LE.  MAXREP)  THEN 

WRITEC 1,3337)  1+70 , INT( DEGFR( 1+70) ) ,CHISQ( 1,1+70), 
+  QHTUPRC 1 , 1+70) ,CHISQC 2 , 1+70) , QHTUPRC  2 , 1+70) 

ELSE 
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-  iW»' 


END  IF  ■*******>-.■ 

i  =  i  +  i 

GOTO  175 
180  END  IF 
ELSE 
ENDIF 

9999  WRITE(1,'("THE  TOTAL  NO  OF  REPS  WAS:",  18)’)  TOTREP 

WRITE(  1, '  (  "THE  TOTAL  NO  OF  EFFECTIVE  REPS  WAS:  "  ,  18)  '  )  LOOP 
WRITE( 1, ' ( ' 'THE  TOTAL  NO  OF  NO  FAILURE  RUNS  WAS: ' ' ,18)' )  ZFAILS 
WRITE(  1 , '  (  '  'THE  TOTAL  NO  OF  RUNS  WITH  FAILURES  WAS:",  18)')  FAILS 
WRITE(  1, '  (  ' '  AVG  NO  OF  COMPONENTS  WITH  NO  FAILURES  PER  REP  WAS:  ", 

+F5. 2) ' )  ZFPREP  /  MAXREP 
0008  FORMAT  (/  3X  ,'C  l'.SX.’C  2’, 

+5X,'C  3' ,5X, 'C  4' ,5X,  C  5',5X,'C  6',5X,’C  7',5X, 

+'C  8’,5X,’C  9',5X,'C  10',4X,’C  ll’,4X, 

+'C  12 ' ,4X, ' C  13',4X,'C  14'  4X,'C  15') 

0009  FORMAT( / IX , ' REP  NO* ,2X. 'DF  ,1X. ’CHISQR( ' ,F4.  3, ' )' ,  IX, 

+  ' QHTUPR( ’ ,F4. 3 , ' ) '  IX, fCHISQR( f  .F4. 3 , ' ) '  . IX, * QHTUPR( ' ,F4. 3 , ' ) ’ , 

+2X , ' REP  NO' , 2X, ' DF  , IX. ' CHISQR(  ,F4. 3 . ' )  , IX, 

+  'QHTUPR(' ,F4.  3,')' ,lX,fCHISQR( ' ,F4.  3, f )' . IX, fQHTUFR( ' ,F4.  3 , ’ ) '/) 

0001  FORMAT  (///'UNIFORM  RANDOM  DEVIATES  ARE:  f) 

0002  FORMAT  (///’BERNOULLI  TRIALS  ARE:  ') 

0003  FORMAT  (///’TOTAL  NO.  OF  FAILURES  FOR  EACH  COMPONENT:  ') 

0004  FORMAT  (///'ESTIMATED  UNRELIABILITY  FOR  EACH  COMPONENT:  ') 

0005  FORMAT  (///'TOTAL  NUMBER  OF  MISSION  TESTS:  ') 

0006  FORMAT  (///'ESTIMATED  WEIGHTS  FOR  EACH  COMPONENT:  ’) 

0007  FORMAT  ( /// ' Q  I  FOR  EACH  COMPONENT: ’) 

1111  FORMAT  (15F8.5) 

2222  FORMAT  (/IX, 15( I4.4X)) 

3333  FORMAT  (/IX, 15(14, 4X)) 

3334  FORMAT  (/15F8.  5) 

3335  FORMAT  ( /IX, 15( 14 ,4X) ) 

3336  FORMAT  (T3 , 14 ,T9 , 13 ,T13 ,F11.  5 ,T27 ,F8.  5 ,T39 ,F11.  5 ,T53 ,F8.  5) 

3337  FORMAT  ( '  +  ' ,T67 , 14 ,T73 , 13 ,T77 ,F11.  5 ,T91 ,F8.  5 ,T103 ,F11.  5 ,T117 ,F8.  5) 

4442  FORMAT  ('  ',///' THE  RESULTING  (1  -  ' ,F4. 3,')  CONFIDENCE  ', 

+  '  LIMIT  IS:  '  ,T50, '  00.  000  '  ) 

4443  FORMAT  ('  ',///' THE  RESULTING  (1  -  ' ,F4. 3,’)  CONFIDENCE  ', 

+  ' LIMIT  IS: ' ,T50,'100.  0000') 

4444  FORMAT  ('  './//'THE  RESULTING  (1  -  ' ,F4. 3,')  CONFIDENCE  ', 

+' LIMIT  IS: ' ,T50 ,F8. 4) 

4445  FORMAT  ('  './’THE  RSHAT  VALUE  CLOSEST  TO  RS  IS:  ’,T51,F8.5) 

4446  FORMAT  ('  './’(FIRST  NEGATIVE  DIFFERENCE)') 

4447  FORMAT  ('  './'(ELEMENT  PRECEEDING  FIRST  NEGATIVE  DIFFERENCE)’) 

4448  FORMAT  ('  './'THE  RHTSTR  VALUE  CLOSEST  TO  RS  IS:  ',T51,F8. 5) 

4449  FORMAT  ('  ’ ./'THE  RSHTBR  VALUE  CLOSEST  TO  RSBRDG  IS:  ',T51,F8. 5) 

5555  FORMAT  ('  ’.///'THE  ' , 14 , ' ( 1- ' ,F4. 3, ' )  QUANTILE  IS:  ' ,T49,F8.  3) 

5556  FORMAT  (’  ’./'THE  VALUE  OF  RSHAT  FOR  THAT  QUANTILE  IS: ',T51,F8.  5) 

5557  FORMAT  (’  ’./'THE  DIFFERENCES  -  RSHAT)  IS:  ’  ,T51  ,F8.  5) 

5656  FORMAT  (’  './'THE  VALUE  OF  RHTSTR  FOR  THAT  QUANTILE  IS:  ' ,T51,F8.5) 

5666  FORMAT  ('  './'THE  VALUE  OF  RSHTBR  FOR  THAT  QUANTILE  IS:  ’,T51,F8.  5) 

5657  FORMAT  ('  ’./'THE  DIFFERENCES  -  RHTSTR)  IS:  '  ,T51  ,F8.  5) 

5667  FORMAT  (’  './'THE  DIFFERENCE( RS  -  RSHTBR)  IS: ' ,T51 ,F6. 5) 

5755  FORMAT  ('  './//'SINCE  THE  NO.  OF  MISSION  TESTS  IS  THE  SAME  FOR', 

+'  ALL  COMPONENTS  THE  CLOSED  FORM  SERIES  SYSTEM  RELIABILITY  ' , 

+  ’ "RHTSTR"  IS  COMPUTED') 


6666  FORMAT  ('+'  1  *********************************************** ' 

+'*  RUN  INPUT  SETTINGS  FOR  ADJUSTED  SUM  OF  THE  NI  AI  ***********' 

4-  *  Vr  Vr  iV  iViVVr  Vr  1V7V  Vr  VctVtVVc  *  ^ 

6667  FORMAT  (’  //'NUMBER  OF  REPLICATIONS: ’ ,T50 , 14) 

6668  FORMAT  ('  //'NUMBER  OF  COMPONENTS:  \TS0,I4) 

6669  FORMAT  ('  ',//' SYSTEM  RELIABILITY  FUNCTION: ’ ,T50, ' SERIES' ) 

6699  FORMAT  ('  //'SYSTEM  RELIABILITY  FUNCTION: ' ,T50 , ' BRIDGE ' ) 

6670  FORMAT  ('  ',//' MASTER  UNRELIABILITY  USED: ' ,T50 ,Fe.  5) 

6770  FORMAT  ('  './/'ZERO  FAILURE  SCALING  FACTOR: ’ .T50.F5.  2) 

6671  FORMAT  ('  ',//'  INPUT  VEIGHTS(A  SUB  l"S):') 

6674  FORMAT  ( '+'  , /// 1 *********************************************** 1 
+  '**r  UN  RESULT  S****************************************' 

+ 1  ************** 1  ) 

5675  FORMAT  (  '+'  f /// 1  ************************************************ 
+  '  ESTIMATE  ERRORS  ********************************************* 1 
***************** ^ 

6676  FORMAT  ( '+'  /// ' ************************************************ 
+'  TRUE  CONFIDENCE  LIMITS  *************************************** 

+ 1  ************** 1 ) 

END 


APPENDIX  C.  INTERVAL  ESTIMATING  PROCEDURE  3 


PROGRAM  ZFYSCY 

^***iV-Vyc***Vr***Vf*****Vf*Vr*yr**yr*iV***Tir****VrTVyryf*yrVr*******iicVf**-**yrT,tyfVir*yr*'ilnlryfVt** 
*  * 

*  TITLE:  BINOMIAL  INTERVAL  ESTIMATION  PROCEDURE  * 

*  ZERO  FAILURES  DISALLOWED;  WITH  SCALING  * 

*  AUTHOR:  E.  F.  BELLINI,  LT,  USN  * 

*  DATE:  NOV  89  * 

*  * 

*  THIS  PROGRAM  COMPUTES  THE  TRUE  CONFIDENCE  LEVEL  FOR  THE  ESTIMATE  * 

*  RELIABILITY  OF  A  SERIES  AND  BRIDGE  SYSTEM  GIVEN  THE  RELIABILITY  * 

*  OF  THEIR  COMPONENTS  * 

*  * 

*  IN  ITS  PRESENT  CONFIGURATION  THIS  PROGRAM  IS  SET  UP  TO  RUN  12  * 

*  TIMES  EACH  TIME  PRODUCING  1000  REPLICATIONS  USING  A  DIFFERENT  * 

*  SET  OF  INPUT  DATA.  RUN  THE  PROGRAM  FROM  CMS  BY  TYPING  1 B 1  EXEC’.* 

*  THE  REXX  EXEC  PROGRAM  * 

*  ’ B 1 ’  CALLS  THE  INPUT  FILES  TO  BE  READ  AND  NAMES  THE  12  OUTPUT  * 

*  FILES  RESULTING  FROM  THE  12  CONSECUTIVE  RUNS.  BY  EDITING  THE 

*  INDEX  COUNTERS  I,  J,  K  OF  THE  ’ B 1 1  EXEC  ONE  CAN  RUN  ANY  USER- 

*  SPECIFIC  RUN  FROM  JUST  ONE  RUN  TO  ALL  12. 

if  if 

*  VARIABLES  USED 


AiiATI 

AI 

ALFA 

B1GF 

CHI5Q 

COUNT 1 

C1C15 

DEGFR 

DELBRG 

DELSTK 

DELTAR 

DIFF 

EPS 

ERROR 

FAILS 

FI 

FLAG 

INC 

KEY1 

KEY  2 

KEY3 

KEY4 

KK 

LOOP 
MAXALF 
MAXREP 
MAX RUN 
MSTKQ 


WEIGHT  ESTIMATES  FOR  EACH  COMPONENT  * 

INPUT  WEIGHTS  FOR  EACH  COMPONENT  * 

LEVELS  OF  SIGNIFICANCE  * 

TOTAL  NO.  OF  FAILURES  FOR  EACH  REPLICATION  * 

CHI-SQUARE  RANDOM  VARIABLE  VALUE  * 

COUNTS  THE  NO.  OF  COMPONENTS  WITH  FAILURES  * 

FORMAT  LABEL  * 

DEGREES  OF  FREEDOM 
DIFFERENCE  FOR  BRIDGE  SYSTEM 

DIFFERENCE  FOR  SERIES  SYSTEM-  CLOSED  FORM  * 

DIFFERENCE  FOR  SERIES  SYSTEM  * 

DIFFERENCE  (TRUE  REL.  -  ESTIMATED  DEL.  ) 

SMALL  QUANTITY! CONSTANT) 

PARAMETER  FOR  IMSL  ROUTINE  * 

COUNTS  NO  REPLICATIONS  WITH  AI  LST.  1  FAILURE  * 
NO.  OF  FAILURES  FOR  EACH  COMPONENT! ALL  MISSION  TST)* 
1  JF  ALL  Cl  -i\  HAVE  SAME  NO.  OF  MISSION  TESTS  * 

INCREMENT  STEP  SIZE  FOR  ROUTINE  USMNMX  * 

ARRAY  OF  INDEXES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT  * 

ARRAY  OF  INDECES  FOR  ROUTINE  SHSORT 
ARRAY  SIZING  PARAMETER  FOR  THE  MAX  NO  OF  COMPONENTS* 
COUNTS  NO.  OF  REPLICATION  PERFORMED 
MAX  NO.  OF  SIGNIFICANT  LEVELS  DESIRED! ARRAY  SIZING)* 
MAX  NO.  OF  REPLICATIONS  * 

MAX  NO.  OF  PROGRAM  ITERATIONS  ALLOWED  * 

MASTER  UNRELIABILITY (USED  WITH  Al’S  TO  CnLC.  Ql’S)  * 


it 

MULT 

* 

N 

it 

NT  MAX 

it 

NIMIM 

* 

NINDX 

i k 

NIREAL 

it 

NMAX 

it 

NPRNT 

it 

PRNT 

it 

QHATI 

it 

QHTMAX 

it 

QHTUPR 

it 

QI 

it 

QINDX 

it 

QUANTL 

it 

REPSHD 

it 

RHTSTR 

it 

RS 

it 

RSBRDG 

it 

RSHAT 

it 

RSHTBR 

it 

SEED 

it 

SELCTA 

it 

SELCTB 

it 

SORT 

it 

SUMNAI 

it 

TEMP 

it 

TOTREP 

* 

TRAN SR 

it 

TRANSQ 

it 

TRANSR 

it 

TRIALS 

it 

TRNSTR 

it 

TRUQNT 

it 

UNIRV 

it 

ZFAILS 

it 

ZFPREP 

it 

MULTIPLIER  FOR  RANDOM  NO.  GENERATOR  SRND 

NO.  OF  MISSION  TEST  FOR  EACH  COMPONENT  * 

MAX  NO.  OF  MISSION  TESTS 

MIN  NO.  OF  MISSION  TESTS 

INDEX  NO.  OF  MAX  NO.  OF  MISSION  TESTS 

NO.  OF  MISSION  TESTS  TRANSFORMED  TO  REAL 

MAX  NO.  OF  MISSION  TESTS  FOR  OUTPUT  CONTROL 

FLAG  FOR  DETAILED  REPORT  OUTPUT 

SAME  AS  ABOVE( PARAMETER) 

UNRELIABILITY  ESTIMATES  FOR  EACH  COMPONENT 
LARGEST  QHATI 

UPPER  LIMIT  ON  SYSTEM  UNRELIABILITY 
INPUT  UNRELIABILIY  FOR  EACH  COMPONENT 
INDEX 
QUANTILE 

REPLICATIONS  HEADING  FORMAT  NUMBER 

SERIES  SYSTEM  RELIABILITY  ESTIMATE( CLOSED  FORM) 

TRUE  SERIES  SYSTEM  RELIABILITY 

TRUE  BRIDGE  SYSTEM  RELIABILITY 

SERIES  SYSTEM  RELIABILITY  ESTIMATE 

BRIDGE  SYSTEM  RELIABILITY  ESTIMATE 

PARAMETER 

SIGNIFICANCE  LEVEL  SELECTION 
SIGNIFICANCE  LEVEL  SELECTION 
PARAMETER  FOR  ROUTINE  SRND 
SUM  OF  THE  PRODUCT  OF  NT'S  AND  Al'S 
TEMPORARY  ARRAY 

TOTAL  NUMBER  OF  PROGRAM  ITERATIONS 

TEMPORARY  ARRAY 

TEMPORARY  ARRAY 

TEMPORARY  ARRAY 

BERNOULLI  TRIALS  ARRAY  (2-DIM) 

TEMPORARY  ARRAY 
TRUE  QUANTILE 

UNIFORM  RANDOM  DEVIATES  (2-DIM) 

TOTAL  NUMBER  OF  REPLICATIONS  WITH  ZERO  FAILURES 
NO.  OF  COMPNTS.  WITH  ZERO  FAILURES  PER  REPLICATION 


* 


* 

* 
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* 

* 
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PARAMETER  (KK=10 ,MAXALF=2 ,NPRNT=0) 

PARAMETER  (MAXREP=1000 ,  MAXRUN=10000 ,  EPS=.  000001) 

REAL*4  UNIRV(  15 , ICOO) ,TEMP( 1000) ,QI(KK) ,  AI(KK),  AHATI(KK) 

REAL*4  QHATI(KK) ,  NMAX,  NNMAX,  QHTMAX ,  ALFA(MAXALF) 

REAL*4  DF(  5 ) , A ALFA ( 5 ) , SUMNAI , RSHAT( MAXALF .MAXREP) , RS 
REAL* 4  KEYl(MAXREP) ,KEY2(MAXREP) ,KEY3( MAXREP) ,TRNSTR( MAXREP ) 
REAL*4  DEGFR( MAXREP),  QHTUPR( MAXALF, MAXREP ) ,CHISQ(MAXALF, MAXREP) 
REAL**  RHTSTR( MAXALF, MAXREP) 

REAL*4  DELTARC MAXALF) ,  TRANSQ( MAXREP) , TRANS R( MAXREP) ,DIFF(MAXRE?) 
REAL* 4  DELSTR( MAXALF),  NN(KK) 

REAL*4  RSHTBR( MAXALF, MAXREP) , DELBRG( MAXALF) , KEY4( MAXREP) 

REAL*4  TRANBR( MAXREP),  RSBRDG 

REAL*4  SCALEN,  MINQI,  MAXQI ,  ZFPREP.MSTRQ 


INTEGER  SEED,  MULT,  SORT,  TRIALS( 15 , 1000) ,  BIGF,  FI(KK) ,  N( KK) 


INTEGER  NINDX,  QINDX,  ERROR,  REPS,  SELCTA,  SELCTQ,  TOTREP 
INTEGER  COUNT 1 ,  C1C15,  REPSHD,  SELCTB ,  ALF,  FLAG,  LOOP.PRNT 
INTEGER  QUANTL(MAXALF) ,  TRUQNT(MAXALF) ,  ZFAILS,  FAILS 
INTEGER  INC 


DATA  SEED/123457/,  MULT/1/,  INC/1/ 

DATA  AALFA/. 01,.  05,.  9,.  95,.  99/,  DF/1, 5, 10,30,40/ 

DATA  ALFA/. 20,. 050/ 

*  DATA  MSTRQ/O. 1,  .05/ 

*  DATA  AI/.  20,.  10,.  05,.  05,.  05/ 

*  DATA  AI/.20,.  10, .05, .05, .05, .05, .05, .05, .05, .05/ 

*  DATA  AI/.  05, .05, .05, .05, .05, .05, .05,. 05, .05, .05/ 

*  DATA  AI/.  05, .05, .05, .05, .05, .05, .05, .05, .05, .05, .05, .05, .05, 

*  +.05,.  05/ 

*  DATA  AI/.  1,.  1 ,.  1,.  1 , 1,.  05 ,.  05,.  2,.  25 ,.  5/ 

*  DATA  AI/.  1,.  1,1.  ,.05,.  2/ 

*  DATA  N/100,100,100,100,100/ 

*  DATA  N/ 10, 6, 5, 2, 10, 4, 5, 6, 7, 2/ 

*  DATA  N/50, 50, 50, 50, 50, 50, 50, 50, 50, 50/ 

*  DATA  N/ 100 , 100 , 100 , 100 , 100 , 100 , 100 , 100 , 100 , 100/ 

*  DATA  N/300, 300, 300, 300, 300, 300, 300, 300, 300, 300/ 

*  DATA  N/20, 20, 20, 20, 20, 20, 20, 20, 20, 20/ 

*  DATA  N/ 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 

*  +1000,1000,1000,1000/ 

DATA  SORT/O/ 

ASSIGN  8  TO  C1C15 
ASSIGN  9  TO  REF 3 HD 


*  CALL  EXCMS( ’FILEDEF  01  DISK  B  OUTPUT  AI  (LRECL  132  ’) 

*  CALL  EXCMSC FILEDEF  02  DISK  JNK  OUTPUT  AI  (LRECL  132  ') 

SELCTQ  =  1 
PRNT  =  NPRNT 

DO  12  1=1, KK 

AI(I)  =  9999. 

N( I)  =  99999999 
12  CONTINUE 

READ(03,*)K,MSTRQ 

DO  11  1=1, K 

READ( 03 ,*)  AI(I),N(I) 

11  CONTINUE 

***/(  INITIALIZE  THE  QHTUPR  ARRAY  OF  UNRELIABILITY  REPLICATIONS,  //*** 

*  RSHAT  ARRAY  OF  ESTIMATE  SYSTEM  RELIABILITY  REPLICATIONS  * 

*  AND  RHTSTR  ARRAY  OF  EST.  SYST.  REL.  FOR  A  SERIES  SYST  WHEN  * 

***//  ALL  THE  COMPONENT  MISSION  TESTS  ARE  EQUAL  IN  NUMBER  //*** 


DO  172  ALF=1 .MAXALF 

DO  173  REPS=1 .MAXREP 

QHTUPRC ALF, REPS)  =  0. 
RSHATC ALF, REPS)  =  0. 
RHTSTRC ALF, REPS)  =  0. 
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RSHTBR(ALF,REPS)  =  0. 

173  CONTINUE 

172  CONTINUE 

***//  SET  FLAG  TO  1  IF  ALL  COMPONENTS  HAVE  SAME  NO.  OF  MISSION  TESTS**** 
FLAG=1 

DO  50  1=1, K  -1 

IF( (N(I)  -  N(I+1)).NE. 0)  THEN 
FLAG=0 

ELSE 
END  IF 
50  CONTINUE 

PRINT  *,  'FLAG  IS: ' ,  FLAG 

***//  MAIN  PROGRAM  OUTER  LOOP  START( EVERY  LOOP  IS  ONE  REPLICATION)//*** 


ZFPREP  =  0. 

FAILS  =  0 
ZFAILS  =  0 
COUNT 1  =  0 
TOTREP  =  0 
LOOP  =  0 

10  IF(LOOP.  LT.MAXREP)  THEN 

IF(TOTREP.  LT.MAXRUN)  THEN 
SELCTA  =  1 
SELCTB  =  2 


ti."V 

***//  FILL  ARRAY  KEY(REPS)  WITH  INTEGERS  1  TO  K  TO  BE  USED  AS  OUTPUT 
***//  OF  THE  SUBROUTINE  SHSORT 

DO  95  REPS=1 ,  MAXREP 
KEYl(REPS)  =  REPS 
KEY2(REPS)  =  REPS 
KEY3(REPS)  =  REPS 
KEY4(REPS)  =  REPS 
95  CONTINUE 

***//  CALCULATE  NMAX  NOT  TO  PRINT  LONGER  THAN  THE  MAX  SAMPLE  SIZE 
***//  CALCULATE  THE  MAXIMUM  NO.  OF  TRIALS  AND  ITS  INDEX  NO.  //*** 

CALL  IMAX( N , K , NMAX , NINDX ) 


***//  CALCULATE  THE  Ql’S  FROM  THE  GIVEN  MASTER  Q  AND  THE  Al'S  //*** 

DO  115  1=1,  K 

QIC  I)  =  MSTRQ  *  AI(I) 

115  CONTINUE 

DO  120  1=1,15 

DO  125  J=l,500 

UNIRV( I , J)  =  999. 

TRIALSC I , J)  =  99999 
125  CONTINUE 

120  CONTINUE 


58 


//*** 


***//  CALCULATE  THE  SCALING  FACTOR  'SCALES’ 

CALL  USMNMX(QI ,K, INC ,MINQI ,MAXQI) 

SCALEN  =  MAXQI  /  MINQI 
***//  DRAW  UNIFORM  (0,1)  RV'S  AND  CONVERT  TO  BERNOULLI  TRIALS  //*** 

DO  130  1=1,  K 

'  CALL  SRND(SEED,  TEMP,  N(I),  MULT,  SORT) 

DO  135  J=l,  N( I) 

UNIRV(T.J)  =  TEMP(J) 

IF  (UNIRV( I , J) . LE.  1  -  QIC  I) )  THEN 
TRIALS( I , J)  =  0 

ELSE 

TRIALS( I , J)  =  1 
END  IF 

135  CONTINUE 
130  CONTINUE 


***//  CALCULATE  THE  NO.  OF  FAILURES  FOR  EACH  COMPONENT  / /*** 

DO  150  1=1,  K 
FI( I)  =  0 
150  CONTINUE 

***//  CALCULATE  THE  F  SUB  I' S  AND  THE  GRAND  TOTAL  NO.  OF  FAILURES 

BIGF  =  0 
DO  155  1=1,  K 

DO  160  J=l,  N( I) 

FI( I)  =  FI C I)  +  TRIALSC I , J) 

160  CONTINUE 

***//  CALCULATE  THE  QHAT  SUB  l'S:  F  SUB  I’S  DIVIDED  BY  N  SUB  I’S 
IF(FI(I). EQ.  0)  THEN 

QHATI(I)  =  1.  /  (SCALEN  *  N(I)) 

ZFPREP  =  ZFPREP  +1 

ELSE 

QHATI(I)  =  REAL(FKI))  /  N(I) 

END  IF 

BIGF  =  BIGF  +  FI( I) 

155  CONTINUE 

***//  FIND  THE  MAX  OF  THE  INDIVIDUAL  COMPONENT  UNRELIABILITIES 
CALL  RMAX(OHATI ,  K,  QHTMAX ,  QINDX) 

IF(PRNT.  EQ.  1)  THEN 
WRITE( 1,0007) 

WRITE( 1 ,C1C15 ) 

WRITE( 1,3334)  QI 
WRITE( 1,0001) 

WRITE( 1 ,C1C15) 

DO  140  J=1 ,NMAX 

WRITE( 1,1111)  (UNIRV(I.J),  1=1,  K) 

140  CONTINUE 
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WRITEC 1,0002) 

WRITE( 1 ,C1C15) 

DO  145  J=1 ,NMAX 

WRITE( 1,2222)  (TRIALS( I , J) ,  1=1,  K) 

145  CONTINUE 

WRITE( 1,0003) 

WRITE( 1 ,C1C15) 

WRITE( 1 , 3333)  FI 
WRITEC 1,0005) 

WRITE( 1 ,C1C15) 

WRITE( 1 , 3335)  N 
WRITEC 1,0004) 

WRITEC 1,C 1C 15) 

-WRITEC 1.33340  QHATi — 

WRITEC 1,’C/  THE  MAXIMUM  Q  HAT  SUB  I  IS:”,  T40.  F8.  5)')  QHTMAX 
WRITEC1,'C/"THE  MAXI  Q  HAT  SUB  I  IS  ELMNT  NO.  :  f '  .T40.I5) '  )  QINDX 
WRITEC  1,’C  /"THE  GRAND  TOTAL  NO.  OF  FAILURES  IS:,f,T40,  15)')  BIGF 
ELSE 
ENDIF 

***//  TEST  FOR  A  REP  WITH  AT  LST  TWO  COMP  WITH  AT  LST  ONE  FAILURE  EACH 
DO  162  J=l,  K 
NN( J)  =  0 

IFCFICJ). GE. 1)  THEN 

COUNT 1  =  COUNT 1  +  1 
NNCJ)=  N( J) 

ELSE 

NNCJ)  =  NCJ)  *  SCALEN 
END  IF 

162  CONTINUE 

IFCCOUNT1.EQ.  0)  THEN 

ZFAILS  =  ZFAILS  +  1 

ELSE 

FAILS  =  FAILS  +  1 
END  IF 

IF( COUNT 1.  GE.  0)  THEN 
LOOP  =  LOOP  +  1 
TOTREP  =  TOTREP  +  1 

ELSE 

TOTREP  =  TOTREP  +  1 
GO  TO  10 
END  IF 

***//  CALCULATE  THE  AHAT  SUB  I'S  C WEIGHT  ESTIMATES) 

SUMNAI  =  0. 

DO  165  1=1,  K 

AKATICI)  =  QHATI(I)  /  QHTMAX 
SUMNAI  =  SUMNAI  +  NNCI)  *  AHATICI) 

165  CONTINUE 

IFCPRNT.EQ.  1)  THEN 

-  WRITEC 1,0006)  — 

WRITEC 1,C1C15) 

WRITEC 1,3334)  AHAT I 

ELSE 
END  IF 
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***//  CALCULATE  1  REPLICATION  OF  UPPR  ALFA  C.  L.  ON  SYSTEM  RELIABILITY 
DEGFR(LOOP)  =  2  *  (1  +  BIGF) 

DO  170  ALF=1 ,  MAXALF 

CALL  MDCHK1  -  ALFA(ALF)  .DEGFR(LOOP)  , CHI SQ( ALF, LOOP)  ,  ERROR) 
QHTUPR( ALF , LOOP)  =  CHISQC ALF .LOOP)  /  (2  *  SUMNAI) 

IFCFLAG.EQ.  1)  THEN 

RHTSTRfALF.LOOPt  =  1  -(CHISCK ALF, LOOP)  /  REAL(2*N(  1))) 

ELSE 
END  IF 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  COMPNTS.  IN  SERIES 

CALL  RHTSRS(QHTUPR(ALF,LOOP),  AHATI.K,  RSHAT( ALF, LOOP) ) 

*  +T40.F8.5)’)  RSHAT( ALF, LOOP) 

***//  CALCULATE  VALUE  OF  THE  SYSTEM  RELIABILITY  FOR  BRIDGE  STRUCTURE  *** 
CALL  RHTBRG(QHTUPR(ALF, LOOP), AHATI ,K,RSHTBR( ALF, LOOP)) 

170  CONTINUE 

***//  THIS  ELSE  AND  ENDIF  ARE  FOR  THE  TEST  AGAINST  MAXRUN  //*** 


WRITE ( 1 , ' (  "  ’ ’  / ’ ' PROGRAM  EXCEEDED  THE  MAX  NO.  OF  RUNS’ ' , 

+  ' ’  ALLOWED  OF:  ’ ’ ,I6)‘)  TDTREP 
GOTO  9999 

END  IF 
GOTO  10 
END  IF 

***//  SORT  THE  ARRAYS  OF  SYSTEM  UNRELIABILITIES( 1  FOR  EACH  CONF.  LEVEL) 


DO  700  ALF=1 ,  MAXALF 

DO  800  REPS=1 .  MAXREP 

TRANSQ(REPS)  =  QHTUPR( ALF, REPS) 
TRANSR(REPS)  =  RSHAT( ALF , REPS) 
TRNSTR(REPS)  =  RHTSTR( ALF, REPS) 
TRANBR(REPS)  =  RSHTBR( ALF, REPS) 
CONTINUE 

CALL  SHS0RT(TRANSQ,KEY1, MAXREP) 

CALL  SHSORT(TRANSR,KEY2, MAXREP) 

CALL  SHSORT(TRNSTR,KEY3, MAXREP) 

CALL  SHS0RT(TRANBR,KEY4, MAXREP) 

~DO  900  REP§=1,  MAXREP 

QHTUPR( ALF, REPS)  =  TRANSQ(REPS) 
RSHAT( ALF, REPS)  =  TRANSR(REPS) 
RHTSTR( ALF, REPS)  =  TRNSTR(REPS) 
RSHTBR( ALF, REPS)  =  TRANBR(REPS) 
CONTINUE 
CONTINUE 
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//  PRINT  OUTPUT  REPORT  HEADINGS 


j  /*** 


WRITEC 1,6666) 

VRITE( 1,6667)  MAXREP 
WRITEC  1,6668)  K  — 

WRITEC 1,6770)  SCALEN 
WRITEC 1,6669) 

IFCK.EQ.  5)  THEN 
WRITEC 1,6699) 

ELSE 
END  IF 

WRITEC 1,6670)  MSTRQ 
WRITEC 1,6671) 

WRITEC 1,C 1C 15) 

WRITEC 1,3334)  AI 
WRITEC 1,0007) 

WRITEC 1,C 1C 15) 

WRITEC 1,3334)  QI 
WRITEC 1,0005) 

WRITEC 1,C 1C 15) 

WRITEC 1,3335)  N 
WRITEC 1,6674) 

WRITEC 2, 6666) 

WRITEC 2, 6667)  MAXREP 
WRITEC 2, 6668)  K 
WRITEC 2, 6770)  SCALEN 
WRITEC 2, 6669) 

IFCK.EQ.  5)  THEN 

WRITEC 1,6699) 

ELSE 
END  IF' 

WRITEC 2, 6670)  MSTRQ 
WRITEC 2,6671) 

WRITEC 2, C1C15) 

WRITEC 2, 3334)  AI 
WRITEC 2, 000 7) 

WRITEC  2 ,C1C15) 

WRITEC2.3334)  QI 
WRITEC 2, 0005) 

WRITEC2 ,C1C15) 

WRITEC 2, 3335)  N 
WRITEC 2, 6674) 

WRITEC  2 , '  C  '  'SORTED  RSHAH  1  IS:  "  ,/L5CF8.  5))  '  ) 
+(RSHAT( 1 , REPS) ,  REPS=1,  MAXREP) 

WRITEC 2,'  C"  SORTED  RSHAT  2  IS:  ’  '  , /15CF8.  5) )  '  ) 

+C RSHATC 2 ,REPS ) ,  REPS=1,  MAXREP) 

IFCFLAG. EQ. 1)  THEN 

WRITE  C  2  ,'C"  SORTED  RHTSTR  1  IS:  "  ,/10CF8.  5))  '  ) 
+C RHTSTRC 1 , REPS ) ,  REPS=1,  MAXREP) 

WRITE C 2, 'C ’ 'SORTED  RHTSTR  2  IS: ' ' ,/10(F8. 5))' ) 
+C RHTSTRC 2, REPS),  REPS=1,  MAXREP) 

ELSE 
END  IF 

IFCK.EQ.  5)  THEN 

WRITEC  2 ,'("  SORTED  RSHTBR  1  IS:  "  ,/10(F8.  5)  )  '  ) 
+CRSHTBRC 1 ,REPS) ,  REPS=1,  MAXREP) 
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-  WRITEC  2 SORTED  RSHT3R  2  ISs-A  ’  ,/10(F8.  5))  '  ) 

+(RSHTBR( 2 ,REPS) ,  REPS=1,  HAXREP)  •  - 

ELSE 
END  IF 

***//  COMPUTE  THE  VALUE  RS  OF  THE  TRUE  SYSTEM  REL.  FNCTN. (SERIES  SYSTEM) 
***//  AND  FOR  THE  5 -COMPONENT  BRIDGE  STRUCTURE 

CALL  RSRS(QI ,K.RS) 

WRITE(  1 , '  (  "  '  '  ,////'  1  THE  TRUE  SERIES  SYSTEM  ”, 

+  ' 'RELIABILITY  VALUE  IS:  ”  ,T51,F8. 5) ’ )  RS 
CALL  RBRIDGCQI ,K, RSBRDG) 

IFCK.EQ.5)  THEN 

WRITE(1,'(”  ”,////” THE  TRUE  BRIDGE  STRUCTURE  ”, 

+  ”  RELIABILITY  VALUE  IS:  ”  ,T51  ,F8.  5)  ’  )  RSBRDG 
ELSE 
END  IF 

WRITEC 1,6675) 

***//  COMPUTE  THE  DIFFERENCE  'DELTAR'  BTWN.  RS  AND  RSHAT  OF  THE  THEO 
***//  RETICAL  QUANTILE  GIVEN  BY  ALFA(MUST .USE  SORTED  RSHAT  ARRAY) 

IF(FLAG.  EQ.  1)  THEN 
WRITE( 1,5755) 

ELSE 

'  END  IF 

DO  450  ALF=1 ,  MAXALF 

QUANTL(ALF)  =  MAXREP  *  (1  -  ALFA(ALF) ) 

DELTAR(ALF)  =  RS  -  RSHAT(ALF,QUANTL( ALF) ) 

DELBRG(ALF)  =  RSBRDG  -  RSHTBR( ALF , QUANTL( ALF) ) 

IF(FLAG. EQ.  1)  THEN 

DELSTR(ALF)  =  RS  -  RHTSTR( ALF , QUANTL( ALF ) ) 

WRITEC 1,5555)  MAXREP,  ALFA (ALF) ,  RE AL( QUANTL( ALF ) ) 
WRITE( 1 ,5656)  RHTSTR( ALF , QUANTL( ALF) ) 

WRITE( 1,5657)  DELSTR(ALF) 

ELSE 
END  IF 

IFCK.EQ.5)  THEN 

DELBRG(ALF)  =-RSBRDG  -  RSHTBRC ALF,QUANTL( ALF) ) 

WRITEC 1,5555)  MAXREP,  ALFA(ALF) ,  REAL(QUANTL( ALF) ) 
WRITEC 1,5666)  RSHTBR( ALF , QUANTLC ALF ) ) 

WRITEC 1,5667)  DELBRG(ALF) 

ELSE 
END  IF 

WRITEC 1,5555)  MAXREP,  ALFA(ALF),  REAL( QUANTLC  ALF) ) 

WRITE ( 1 , 5 5 5  6 )  RSHAT( ALF , QUANTLC  ALF ) ) 

WRITEC 1,5557)  DELTAR(ALF) 

450  CONTINUE 

***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 
*** j  j  ********  RSHAT  *******  j /*** 

WRITEC 1,6676) 

DO  400  ALF=1, MAXALF 
TRUQNTC ALF)  =0 

DO  500  REPS=1 ,  MAXREP 

DIFF(REPS)  =  RS  -  RSHATC ALF, REPS) 
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CONTINUE 

DO  600  REPS=1,  MAXREP 

IF( ABS(DIFF(REPS) ) . LE. EPS)  THEN 
TRUQNT(ALF)  =  REPS 

WRITE( 1 , ' (  "  "  , /  " TRUE  CONFIDENCE  LIMIT  IS: ' ' , 

F8.  4) '  ) 

(TRUQNT(ALF)  /  REAL(MAXREP) )  *  100. 

GO  TO  620 

ELSEIF(DIFF(REPS).  LT.  0.  )  THEN 
TRUQNT(ALF)  =  REPS 
GO  TO  610 

ELSE 
END  IF 
CONTINUE 

IF(TRUQNT(ALF).EQ.  0.  )  THEN 
WRITEC 1,4443)  ALFA(ALF) 

WRITE(  1,  *  (  1 1  ",/'  'THE  SMALLEST’ '  , 

"  DIFFERENCE  EETWEEN  RS  AND  RSHAT  IS:  "  ,F10.  5) '  )  DIFF( 
MAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1.  )  THEN 
WRITE(  1,4442)  ALFA(ALF) 

WRITE ( 1 , ' (  "  "  ,  /  "  ALL  RSHAT' ' , 
h  "ARE  GREATER  THAN  RS  "  ) ' ) 

ELSEIF(ABS(DIFF(TRUQNT(ALF) )). LE. ABS(DIFF(TRUQNT( ALF)  -  1))) 
h  THEN 

WRITE (1,4444)  ALF A( ALF), 
h  (TRUQNT(ALF)  /  REAL( MAXREP ) )  *  100. 

WRITEC 1,4445)  RSHATC ALF,TRUQNT(ALF) ) 

WRITF.C  1,4446) 

ELSE 

WRITEC 1,4444)  ALFA(ALF) , 

( ( TRUQNT ( ALF )  - 1 )  /  REAL( MAXREP )  )  *  100. 

WRITE ( 1 , 4445 )  RSHAT ( ALF , TRUQNT ( ALF) - 1 ) 

WRITEC 1,4447) 

END  IF 
CONTINUE 


***//  FIND  the  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //*** 

***//  ********  RSHTBR  (BRIDGE)  *******  //*** 


IFCK.EQ.  5)  THEN 
DO  401  ALF=1 .MAXALF 
TRUQNT (ALF)  =  0 

DO  501  REPS=1 ,  MAXREP 

DIFF(REPS)  =  RSBRDG  -  RSHTBRC ALF, REPS) 

CONTINUE 

DO  601  REPS=1,  MAXREP 

IF( ABS(DIFF(REPS)).  LE. EPS)  THEN 
TRUQNT (ALF)  =  REPS 

WRITEC  1,'f"  "  ,  /  "  TRUE  CONFIDENCE  LIMIT  IS:' 
«-  F8.  4)  ' ) 

f  (TRUQNT(ALF)  /  REAL( MAXREP) )  *  100. 

GO  TO  621 

ELSEIF(DIFF(REPS).  LT.  0.  )  THEN 
TRUQNT(ALF)  =  REPS 
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GO  TO  611 

ELSE 
END  IF 

601  CONTINUE 

611  IF(TRUQNT(ALF).  EQ.  0.  )  THEN 

WRITEC 1,4443)  ALFA(ALF) 

WRITE(1,'("  "  ./"THE  SMALLEST"  , 

+  ' '  DIFFERENCE  BETWEEN  RSBRDG  AND  RSHTBR  IS:  "  , 

+  F10. 5) ' )  DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).  EQ.  1.  )  THEN 
WRITE (1,4442)  ALFA(ALF) 

WRITE(1,'("  ",/"  ALL  RSHTBR"  , 

+  "ARE  GREATER  THAN  RSBRDG' ' )' ) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF)  -  1))) 
+  THEN 

WRITE( 1,4444)  ALFA(ALF), 

,  +  ’*  ■,  (TRUQNT(ALF)  /  REAL(MAXREP) )  *  100. 

WRITE( 1,4449)  RSHTBR ( ALF , TRUQNT( ALF ) ) 

WRITEC 1,4446) 

ELSE 


WRITEC 1,4444)  ALFA(ALF), 

+  ((TRUQNT(ALF)-l)  /  REAL(MAXREP))  *  100. 

WRITE( 1,4449)  RSHTBR ( ALF ,TRUQNT( ALF) - 1 ) 
WRITEC 1,4447) 

621  END  IF 

401  CONTINUE 
ELSE 
END  IF 


■mi-  •-«** 

■***//  FIND  THE  TRUE  CONFIDENCE  LEVEL  OF  THE  SYSTEM  REL.  ESTIMATE  //***  ' 
***  j  J  ********  RHTSTR  ********  j  /*** 


IFCFLAG. EQ. 1)  THEN 
DO  4400  ALF=1 ,MAXALF 
TRUQNTC ALF)  =  0 

DO  5500  REPS=1,  MAXREP 

DIFF(REPS)  =  RS  -  RHTSTR (ALF, REPS) 

5500  CONTINUE 

DO  6600  REPS=1,  MAXREP 

IF(ABS(DIFF(REPS) ). LE. EPS)  THEN 
TRUQNTC ALF)  =  REPS 

WRITEC  1,'("  "  , / " TRUE  CONFIDENCE  LIMIT  IS: 
+  F8. 4) ' ) 

+  (TRUQNTC ALF)  /  REAL( MAXREP) )  *  100. 

GO  TO  6620 

ELSEIFCDIFF(REPS). LT. 0. )  THEN 
TRUQNTC ALF)  =  REPS 
GO  TO  6610 

ELSE 
END  IF 

6600  CONTINUE 


6610  IFCTRUQNTC ALF). EQ. 0. )  THEN 

WRITEC 1,4443)  ALFA (ALF) 

WRITEC  1,’C"  ' './"THE  SMALLEST",  • 

+  ' '  DIFFERENCE  BETWEEN  RS  AND  RHTSTR  IS:  "  , 


I  I 


> 
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+  F9. 5)')  DIFF(MAXREP) 

ELSEIF(TRUQNT( ALF).  EQ.  1.  )  THEN 
VRITEC 1,4442)  ALFA (ALF) 

WRITECl.'C"  "./'’ALL  RHTSTR", 

+  "ARE  GREATER  THAN  RS  "  ) ' ) 

ELSE IF(ABS(DIFF(TRUQNT( ALF) )). LE. ABS(DIFF(TRUQNT(ALF)  -  1))) 
+  THEN 

WRITE( 1,4444)  ALFA(ALF) , 

+  (TRUQNT(ALF)  /  REAL(MAXREP) )  *  100. 

VRITE( 1,4448)  RHTSTR ( ALF, TRUQNT( ALF)) 

VRITEC 1,4446) 

ELSE 

VRITEC  1,4444)  ALFA(ALF) , 

+  ( (TRUQNTC ALF) -1)  /  REAL(MAXREP))  *  100. 

VRITEC 1 ,4448 )  RHTSTR ( ALF , TRUQNTC ALF) - 1) 

VRITEC 1,4447) 

6620  END  IF 

4400  CONTINUE 
ELSE 
END  IF 


***//  PRINT  THE  ARRAYS  PERTINENT  TO  THE  OUPUT  OF  EACH  REPLICATION  //*** 
IFCPRNT.EQ. 1)  THEN 

185  VRITEC l.REPSHD)  ALFA(SELCTA) ,  ALFA( SELCTA) , 

+ALFA(SELCTB) , ALFA( SELCTB) ,ALFA( SELCTA) ,ALFA( SELCTA) , ALFAf SELCTB) , 
+ALFA(SELCTB) 

175  IFCI.GE. (MAXREP  +  1))  THEN 
GOTO  180 

ELSE 

IF(  (I.  EQ.  71).  OR.  ( I.  EQ.  211).  OR.  (I.  EQ.  351).  OR.  (I.  EQ.  491).  OR. 

+  (I  -  EQ-  631).  OR.  ( I.  EQ.  7  71).  OR.  (I.  EQ.  911).  OR.  (I.  EQ.  1051)  )  THEN 
I  =  I  +  70 
WRITECl.’C'  '  +  ")') 

GOTO  185 

ELSE 

VRITEC 1,3336)  I,  INT(DEGFR( I) ) ,  CHISQ(1,I),  QHTUPR(l.I), 

+  CHISQC 2,1),  QHTUPRC 2,1) 

END  IF 

IF( I  +  70. LE. MAXREP)  THEN 

VRITEC 1 , 3337)  1+70 . INT(DEGFR( 1+70) ) , CHISQC 1 , 1+70) , 

+  QHTUPRC 1 , 1+70) ,CHISQ( 2 , 1+70) , QHTUPRC  2 , 1+70) 

ELSE 
END  IF 
1  =  1+1 
GOTO  175 
180  END  IF 
ELSE' 

END  IF 

9999  VRITEC  1,'(  "THE  TOTAL  NO  OF  REPS  WAS:  ",I8)’)  TOTREP 

VRITEC  1,'(" THE  TOTAL  NO  OF  EFFECTIVE  REPS  WAS:",  18)')  LOOP 
VRITEC  1,'(" THE  TOTAL  NO  OF  NO  FAILURE  RUNS  WAS:  ",I8)’)  ZFAILS 
VRITEC  1,'(  "THE  TOTAL  NO  OF  RUNS  WITH  FAILURES'  WAS:  ",  18)  ’  )  FAILS 
VRITEC  1  ,*(" AVG  NO  OF  COMPONENTS  WITH  NO  FAILURES  PER  REP  WAS:  ", 
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-i  r  5.  2) '  )  2FPREP  /  MAXREP 
0C08  FORMAT  (/  3X  'C  l',5X.'C  2', 

+5X,'C  3',5X,fC  4' ,5X,  C  5',5X,'C  6’,5X,'C  7',5X, 

+  C  8  .5X,  C  9’,5X,'C  10'  ,4X, ' C  ll\4X, 

+'C  12f,4X,'C  13',4X,'C  14', 4X,'C  15') 

0009  FORMAT( /  IX , ' REP  NO' ,2X  ' DFf , IX, ' CHISQR( ' ,F4.  3, ’ ) ' , IX, 

+  QKTUPR( ' ,F4.  3, ' ) '  IX, fCHISQR( f .F4. 3, '  )  ’  IX, 'QHTUPR( { ,F4. 3, ' ) ' , 
+2X,  REP  NO  ,2X,  DF  ,1X. ’CHISQR(  ,F4. 3, '  )  ,1X, 

+  'QHTUPR( ' ,F4.  3, ' )' ,1X, (CHISQR( ' ,F4.  3,  )' ,1X, fQHTUPR( ' ,F4.  3,')’/) 
0001  FORMAT  (///'UNIFORM  RANDOM  DEVIATES  ARE: ') 

0002  FORMAT  (///'BERNOULLI  TRIALS  ARE:’) 

0003  FORMAT  (///'TOTAL  NO.  OF  FAILURES  FOR  EACH  COMPONENT: ') 

0004  FORMAT  (///'ESTIMATED  UNRELIABILITY  FOR  EACH  COMPONENT:') 

0005  FORMAT  (///'TOTAL  NUMBER  OF  MISSION  TESTS: ') 

0006  FORMAT  (///'ESTIMATED  WEIGHTS  FOR  EACH  COMPONENT: ') 

0007  FORMAT  ( /// ' Q  I  FOR  EACH  COMPONENT:') 

1111  FORMAT  (15F8.5) 

2222  FORMAT  ( /IX, 15( I4,4X) ) 

3333  FORMAT”‘(  / IX ,  15 (  I4y4X) )  ^ 

3334  FORMAT  (/15F8.5) 

3335  FORMAT  ( /IX, 15( 14 ,4X) ) 

3336  FORMAT  (T3 , 14 ,T9 , 13 ,T13 ,F11.  5 ,T27 ,F8.  5 ,T39 ,F11. 5 ,T53 ,F8. 5) 

3337  FORMAT  (’  +  ' ,T67 . 14 ,T73 , 13 ,T77 ,F11. 5 .T91 ,F8. 5 ,T103 ,F11.  5 ,T117 ,F8.  5) 

4442  FORMAT  ('  './//'THE  RESULTING  (1  -  '.F4.3,')  CONFIDENCE 
+' LIMIT  IS: ' ,T50  '  00.000  ') 

4443  FORMAT  ('  './//'THE  RESULTING  (1  -  ' ,F4. 3,’)  CONFIDENCE  ', 

+' LIMIT  IS: ' ,T50  ’ 100. 0000’ ) 

4444  FORMAT  ('  ',///fTHE  RESULTING  (1  -  ’ ,F4. 3,')  CONFIDENCE  ’ , 

+'LIMIT  IS: ' ,T50,F8. 4) 

4445  FORMAT  ('  ',/'THE  RSHAT  VALUE  CLOSEST  TO  RS  IS:  ’iT51,F8. 5) 

4446  FORMAT  ('  ', /’(FIRST  NEGATIVE  DIFFERENCE)’) 

4447  FORMAT  ('  './'(ELEMENT  PRECEEDING  FIRST  NEGATIVE  DIFFERENCE)’) 

4448  FORMAT  ('  ' ,/'THE  RHTSTR  VALUE  CLOSEST  TO  RS  IS:  ’,T51,F8.5) 

4449  FORMAT  (’  ’./’THE  RSHTBR  VALUE  CLOSEST  TO  RSBRDG  IS:  ’,T51,F8.5) 

5555  FORMAT  ('  ’,///’ THE  ' ,14, ' ( 1-' ,F4. 3, ' )  QUANTILE  IS: ' ,T49 ,F8. 3) 

5556  FORMAT  ('  './'THE  VALUE  OF  RSHAT  FOR  THAT  QUANTILE  IS: ' ,T51,F8.  5) 

5557  FORMAT  ('  './'THE  DIFFERENCE(RS  -  RSHAT)  IS: ' ,T51 ,F8. 5) 

5656  FORMAT  ('  ',/'THE  VALUE  OF  RHTSTR  FOR  THAT  QUANTILE  IS:  ',T51,F8.  5) 

5666  FORMAT  ('  ' ,/’THE  VALUE  OF  RSHTBR  FOR  THAT  QUANTILE  IS: ',T51,F8. 5) 

5657  FORMAT  ('  './'THE  DIFFERENCE( RS  -  RHTSTR)  IS: ' ,T51 ,F8. 5) 

5667  FORMAT  ('  './’THE  DIFFERENCES  -  RSHTBR)  IS:  '  ,T51  ,F8.  5) 

5755  FORMAT  ('  './//'SINCE  THE  NO.  OF  MISSION  TESTS  IS  THE  SAME  FOR’, 

+’  ALL  COMPONENTS  THE  CLOSED  FORM  SERIES  SYSTEM  RELIABILITY  ' , 

+  '  "RHTSTR'  ;  IS  COMPUTED') 

+'*  RUN  INPUT  SETTINGS  FOR  ADJUSTED  SUM  OF  THE  NI  AI  ***********' , 

4. 1  ir  *  i\-  * -V  i;  V:  i;  i:  ic  ic  *.'r  *  *  *  Vr  Vr  Vr  -,V  i:  1  ^ 

6667  FORMAT  ('  './/'NUMBER  OF  REPLICATIONS: ' ,T50 , 14) 

6668  FORMAT  ('  './/'NUMBER  OF  COMPONENTS: ' ,T5 0,1 4) 

6669  FORMAT  (’  './/'SYSTEM  RELIABILITY  FUNCTION: ’ ,T50 ,' SERIES' ) 

6699  FORMAT  (’  ’, //’SYSTEM  RELIABILITY  FUNCTION: ’ ,T50 ,' BRIDGE ' ) 

6670  FORMAT  (’  './/'MASTER  UNRELIABILITY  USED:  ' .T50.F8.  5) 

6770  FORMAT  ('  ’.//'ZERO  FAILURE  SCALING  FACTOR: ' ,T50 ,F5. 2) 

6671  FORMAT  ('  '  , // '  INPUT  WF.IGHTS(  A  SUB  I  ”  S) :  '  ) 

6674  FORMAT  (  '+'  ,//  /  1  ' 


+' LIMIT  IS: ' 

4443  FORMAT  ('  ' 
+' LIMIT  IS: ' 

4444  FORMAT  ('  ' 
+' LIMIT  IS: ' 

4445  FORMAT  ('  ' 


5656  FORMAT  ('  ' 

5666  FORMAT  ('  ' 

5657  FORMAT  ('  ' 

5667  FORMAT  ('  ' 
5755  FORMAT  (’  ' 
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+  !  UN  RESULT  S*****************************^**********  ’ 

+  1  **  ************  *  ) 

6675  FORMAT  (  1 +  *  y  ///  *  ************************************************ 
+  *  ESTIMATE  ERRORS  **■******************************************!> ' 

4-  1  **************  1  ) 

66  7  6  FORMAT  (  '  +  1  ,  /  /  /  *  *********************************************** ' 
+ 1  TRUE  CONFIDENCE  LIMITS  *4************************************* 

+  1  ************** 1  } 

END 


APPENDIX  D.  EXTERNAL  SUBROUTINES 


These  six  subroutines  are  used  by  all  three  programs  listed  in  appendices  one 

through  three.  They  must  be  appended  to  the  particular  main  program  being  run  or  they 

must  be  available  on  the  same  disk  as  the  one  from  which  the  main  program  is  being  run. 
SUBROUTINE  IMAX(SEQ,  N,  MX,  INDEX) 

***//  THIS  ROUTINE  IDENTIFIES  THE  MAXIMUM  ELEMENT  OF  AN  INTEGER  VECTOR 

*  NOTE  THAT  SINCE  THE  TEST  IS  .  GT.  SUBROUTINE  ONLY  PICKS  THE  FIRST 

*  OCCURRENCE  OF  A  MAX  SUCH  AS  IN  THE  CASE  OF  A  TIE. 

***//  ALSO,  ARRAY  PASSED  MUST  BE  TYPE  INTEGER 

REAL*4  MX 

INTEGER  SEQ(N),  N,  I,  INDEX 
MX  =  0. 

DO  5  1=1,  N 

IF(SEQ( I). GT. MX)  THEN 
MX  =  SEQ(I) 

INDEX  =  I 

ELSE 
END  IF 

5  CONTINUE 
END 


SUBROUTINE  RBRIDG(QI ,N,RRSS) 

***//  THIS  SUBROUTINE  CALCULATES  THE  "TRUE"  RELIABILITY  OF  A  5 -COMPONENT 
***//  BRIDGE  STRUCTURE 

REAL*4  QI(N),  RRSS 
INTEGER  N 
IF(N.NE.S)  THEN 

WRITEC1,' ("WARNING:  BRIDGE  STRUCTURE  ONLY  USES  ", 

+  '  'THE  FIRST  5  COMPONENTS")') 

ELSE 
END  IF 

RRSS=( 1-QI( 1) )*( 1-QI(4) )+( 1-QI( 2) )*( 1-QI(5) )+( 1-QI( 1) )*( 1-QI( 3) )* 
C(l-QI(S))+( 1-QI(2) )*( 1-QI (3) )*( 1-QI(4) )  -( 1-QI( 1) )*( 1-QI( 2) )* 

C( 1-QI( 3) )*( 1-QI(4) )-( 1-QI( 1) )*( 1-QI( 2) )*( 1-QI( 3) )*( 1-QI(5) )- 
C( 1-QI( 1) )*( 1-QI( 2) )*( 1-QI (4) )*(l-QI(5))-( 1-QI ( 1) )*( 1-QI( 3))* 

C( 1-QI(4) )*( 1-QI(5) ) -( 1-QI( 2) )*( 1-QI( 3))*( 1-QI(4) )*( 1-QI(5))+ 

C2*( 1-QI( 1) )*( 1-QI( 2))*( 1-QI(3))*( 1-QI(4) )*( 1-QI(5)) 

END 


SUBROUTINE  RHTBRG( QHTUP , AHT , N , RRBRDG ) 

***//  THIS  SUBROUTINE  CALCULATES  THE  ESTIMATED  RELIABILITY  OF  A 
***//  5-COMPONENT  BRIDGE  STRUCTURE.  (ONLY  CARRIED  OUT  TO  THE  Q-CUBED  TERM 
REAL*4  QHTUP,  RRBRDG,  AHT(N) 

INTEGER  N 

RRBRDG=  1  -  (  (  QHTUP**  2  )  *(  AHT(  1 )  *  AHT(  2  )  +AHT( 4)*AHT(  5  )  )  )  - 
C( ( QHTUP** 3 ) *( AHT ( 1)*A1IT( 3)*AHT( 5 )+AHT( 2)*AHT( 3)*AHT( 4) ) )+ 
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C ( ( QKTUP,V,V4  )*(  AKT(  1 ) *AHT(  2 )  *AHTC  3 ) *AHT(  4  )  +AHT(  1  )*AHT(  2  ) *AHT(  3 )  * 
CAHT(  5 )+AHT( 1 )*AHT( 2 )*AHT( 4)*AHT( 5 )+AHT( 1 )*AHT( 3 )*AHT( 4)*AHT( 5 )  + 
CAHT( 2 )*AHT( 3 )*AHT( 4)*AHT( 5 ) ) )  - 

C2* ( ( QHTUP**5  )*(AHT(  l ) *AHT(  2 )  *AHT(  3 ) *AHT(  4 ) *AHT(  5 ) ) ) 

END 


SUBROUTINE  RHTSRS(QHTUP, AAHTI ,N,RRSHAT) 

***//  THIS  ROUTINE  CALCULATES  THE  VALUE  OF  THE  SYSTEM  RELIABILITY  OF  A 
*  SERIES  SYSTEM  OF  ' N '  NO.  OF  COMPONENTS  WHICH  HAVE  UNRELIABILITY 
***//  ' QKTUP ' .  THE  FINAL  SYSTEM  RELIABILITY  VALUE  PASSED  IS  ' RRSHAT' 
R£AL*4  QHTUP,  RRSHAT,  AAHTI(N) 

INTEGER  I,  N 
RRSHAT  =  1. 

DO  100  1=1,  N 

RRSHAT  =  RRSHAT  *  Cl  -(AAHTI (I)*  QHTUP) ) 

100  CONTINUE 
END 


SUBROUTINE  RMAX(SEQ,  N,  MX,  INDEX) 

***//  THID  ROUTINE  IDENTIFIES  THE  MAXIMUM  ELEMENT  OF  A  REAL  VECTOR 

*  NOTE  THAT  SINCE  THE  TEST  IS  . GT.  SUBROUTINE  ONLY  PICKS  THE  FIRST 

*  OCCURRENCE  OF  A  MAX  SUCH  AG  IN  THE  CASE  OF  A  TIE. 

***//  ALSO,  ARRAY  PASSED  MUST  BE  TVPE  REAL 

REAL*4  MX,  SEQ(N) 

INTEGER  N,  I,  INDEX 
MX  =  0. 

DO  5  1=1,  N 

IF(SEQ( I).  GT. MX)  THEN 
MX  =  SEQ(I) 

INDEX  =  I 

ELSE 
END  IF 

5  CONTINUE 
END 


SUBROUTINE  RSRS(QIS,N,RRS) 

***//  THIS  ROUTINE  CALCULATES  THE  VALUE  OF  THE  SYSTEM  RELIABILITY  OF  A 
*  SERIES  SYSTEM  OF  'n'  COMPONENTS  WHICH  HAVE  UNRELIABILITY 
***//  1 QIS ' .  THE  FINAL  SYSTEM  RELIABILITY  VALUE  PASSED  IS  'RRS' 

RF,AL*4  QIS(N),  RRS 
INTEGER  I,  N 
RRS  =  1. 

DO  100  1=1,  N 

RRS  =  RRS  *  (1  -  QIS(I)) 

100  CONTINUE 
END 


T'J 
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